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1.0  INTRODUCTION 


As  radar  systems  designs  become  more  complex  and  versatile,  their  per¬ 
formance  becomes  increasingly  sensitive  to  the  operational  environment.  This 
places  an  Increased  burden  on  the  designer  to  Incorporate  the  effects  of  the 
envltonment  In  system  design  studies.  However,  before  this  can  be  done,  it  Is 
necessary  to  develop  an  accurate  model  of  the  environment.  In  the  case  of 
ground  clutter  or  multipath,  this  means  that  there  is  a  need  for  rough  sur¬ 
face  scattering  models.  Such  models  must  nut  only  be  based  upon  sound  physi¬ 
cal  principles  but  also  exhibit  agreement  with  measurements. 

The  purpose  of  this  study  la  to  provide  improved  models  for  surface  scat¬ 
tering  from  terrain  in  the  microwave  frequency  range  and  near  grazing  inci¬ 
dence.  The  basic  approach  entails  applying  the  composite  surface  scattering 
theory  to  a  lossy,  rough  dielectric  surface  .  As  long  as  the  Inci¬ 
dent  angle  is  not  too  near  grazing  and  the  surface  is  reasonably  free  of  sharp 
edges  or  cusps,  the  composite  model  should  be  a  reasonable  description  of  the 
scattering  process.  Since  the  composite  model  Is  based  upon  the  combining  of 
two  asymptotic  scattering  theories.  It  is  approximate.  Thus,  an  additional 
goal  of  this  study  is  to  investigate  new  techniques  for  improving  the  compos¬ 
ite  model.  As  a  first  step  toward  obtaining  improvements  to  the  composite 
model,  a  rigorous  new  formulation  of  the  problem  of  coherent  scattering  from 
a  rough  surface  Is  developed.  The  intent  of  this  work,  is  an  attempt  to  gain 
a  better  understanding  of  the  Interplay  between  the  statistical  surface  param¬ 
eters  and  the  scattered  field.  Such  understanding  is  an  absolutely  essential 
prerequisite  to  modeling  more  complex  factors  such  as  vegetation  and  snow 


cover. 


I.l  Summary  of  Reaulta 


Section  2  corrects  an  error  in  th*’  composite  surface,  scattering  model. 

In  particular,  It  Is  shovm  that  shadowing  is  improperly  accounted  for  and  this 
error  is  corrected.  In  the  corrected  version  of  the  composite  model,  it  is 
shown  that  the  -onventlonal  shadowing  function  multiplies  both  the  zeroth  and 
first  order  incoherent  scattered  power  perturbation  terms.  Thus,  even  for  a 
perfectly  conducting  surface  the  first  order  term  will  go  to  zero  for  back- 
scattering  at  grazing  incidence  due  to  the  shadowing  function.  The  first 
order  perturbation  power  suffers  an  additional  attenuation  due  to  the  shadow¬ 
ing  of  unfavorably  oriented  large  scale  surface  slopes;  however,  this  effect 
is  relatively  small  compared  to  the  impact  of  the  conventional  shadowing  func¬ 
tion. 

Section  3  demonstrates  how  the  shadowing  function  for  non-Gaussian  surfaces 
may  he  easily  obtained  from  existing  shadowing  theories.  The  important 
surface  characteristic  in  the  general  case  is  the  probability  density  func¬ 
tion  of  the  large  scale  slopes  in  the  plane  of  incidence*  Explicit  results 
are  obtained  for  a  surface  characterized  by  a  roughness  whose  probability  den¬ 
sity  function  is  exponential.  Tine  results  of  this  study  are  particularly  Im¬ 
portant  for  terrain  scattering  because  terrain  height  cannot  always  be  described 
by  a  Gaussian  probability  density  function. 

Section  4  extends  the  composite  model  to  bistatic  scattering  from  a 
lossy,  dielectric,  rough  surface.  The  details  are  presented  for  three  cases 
of  Increasing  complexity;  backscatterlng  from  a  surface  with  only  small  scale 
roughness,  bistatlc  scattering  from  a  surface  with  only  small  scale  roughness, 
and,  finally,  bistatic  scattering  from  a  composite  (large  and  small  scales  of 
roughness)  surface.  This  approach  is  a  logical  progression  from  the  simple 
to  the  complex  and  is  therefore  beneficial  to  the  reader.  Furthermore,  this 
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approach  facilitates  checking  the  results  of  the  perturbation  theory  against 
existing  solutions. 

Section  5  discusses  one  technique  for  relating  the  Joint  probability 
density  function  of  the  surface  heights  to  the  joint  density  function  for  the 
slopes.  The  technique  was  originally  obtained  from  an  analysis  of  optical 
scattering  from  a  rough  surface  but  its  relevance  to  this  problem  has  apparently 
been  overlooked.  Although  the  technique  is  not  always  applicable  to  measured 
data,  there  are  cases  where  It  can  provide  the  desired  transformation. 

Section  6  develops  a  new  approach  to  the  problem  of  coherent  scattering 
from  a  perfectly  conducting  rough  surface  based  upon  the  magnetic  field  inte¬ 
gral  equation.  In  contrast  with  the  classical  multiple  scattering  formalism 
which  leads  to  an  Infinite  number  of  integral  equations,  this  approach  results 
in  a  single  integral  equation  of  infinite  dimension.  The  infinite  dimension¬ 
ality  Is  a  consequence  of  retaining  all  orders  of  surface  height  derivatives 
in  the  averaging  process.  The  major  benefit  of  this  approach  is  that  it  is 
possible  to  put  the  mathematical  operation  of  truncating  the  dimensionality 
of  the  Integral  equation  into  one-to-one  correspondence  with  the  neglect  of 
higher  order  surface  height  derivatives.  Comparisons  with  the  multiple  scat¬ 
tering  approach  results  show  very  good  agreement  in  domains  where  both  theories 
are  valid.  This  approach  also  shows  that  in  order  to  neglect  surface  slopes 
the  product  of  the  Rayleigh  roughness  parameter  and  the  rms  surface  slope  must 
be  much  less  than  unity;  that  is  it  is  not  sufficient  to  simply  require  the 
mean  square  slope  to  be  small. 
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2.0  A  CORRECTION  TO  THE  COMPOSITE  SURFACE  SCATTERING  MODEL 


2.1  Background 

In  [1],  a  solution  to  the  problem  of  backscatterlng  from  a  randomly  rough, 
perfectly  conducting  surface  comprising  both  large  and  si  11  scales  of  rough  • 
ness  was  presented.  As  a  direct  consequence  of  the  stlpuiaulon  that  the  sur¬ 
face  roughness  was  a  zero  mean  jointly  Gaussian  process,  the  scattering  cross 
section  per  unit  area  was  determined  to  be  the  sum  of  two  terms,  l.e. 

Op°p»(04)  “  t%p*  is 

dominant  near  normal  Incidence  and  results  from  the  shadow  corrected  optical 
like  reflection  from  properly  oriented  facets  or  specular  points  on  the  sur¬ 
face.  The  [Opp,  (Q,*}*)  term  is  due  to  Bragg  resonance  scattering  from  the 
small  scale  surface  features  with  appropriate  accounting  for  the  resonance 
broadening  e''fects  of  the  large  scale  surface  undulations.  Although  shadowing 
is  formally  accounted  for  in  a  correct  manner  in  the 

there  is  an  error  in  the  exact  representation  of  the  shadowing  function  which 
leads  to  an  incorrect  estimate  of  the  effects  of  shadowing  on  large  angle  of 
incidence  scattering.  The  goal  of  this  section  is  to  correct  the  above  error 
and  to  properly  account  for  the  effect  of  large  scale  shadowing  on  the  small 
scale  scattering  term. 


2.2  Discussion  of  the  Error 

The  analysis  presented  in  [1]  relating  to  the  determination  of  tcfpp,  (G,({))  ] 
is  correct  up  to  and  including  equation  (23).  The  problem  with  the  analysis 
following  (23)  is  a  result  of  inadequate  attention  to  the  definition  of  the 
shadowing  function  R(D,<|))  .  That  is,  R(0,(j))  ,  as  it  appears  in  (23)  of  [1], 
can  and  should  he  expressed  as  follows; 


R(0.(J)) 


,(o) 


I  "-Hx 


Jix  r 

B  ’ 


k  \ 

-f) 


where 


^(^1  I  I  *  ^Ay)X 

H 

la  fhe  probability  that  an  Incident  ray  having  direction  (t^  will  Intersect  a 

point  on  the  large  scale  surface  with  orthogonal  slopes  Cj^j^  and  arid  will 

not  be  shadowed  by  any  other  part  of  the  surface  regardless  of  the  height 

of  the  point  in  question.  The  symbol  in  (2.2)  denotes  the  ensemble 

^  /  \ 

average  over  all  values  of  large  scale  height  and  the  notation  P'^  '  {*) 

Is  the  same  as  that  employed  by  Sancer  (21  in  his  excellent  analysis  of  the 
effect  of  shadowing  on  l^pp»  Iq  '  should  be  noted  that  as  a  result 
of  the  Integrations  in  (23)  of  [1],  the  shadowing  function  R(6,<1))  is  equal 
to  P^°^(*)  evaluated  at  the  specific  large  scale  slope  values  given  by 
Co  ■  k  /B  and  Ca  ■  k  /B  where  k  “  -2k  sin  6  cos  i}>  »  k  ■  -2k  sin  9  sin  <)>  , 
B  ■  2k^co8  0,9  is  the  angle  of  Incidence  relative  to  the  normal  to  the  mean 
(^•0)  plane,  and  4i  is  the  azimuth  direction  of  incidence. 

In  order  to  more  clearly  understand  and  therefore  rectify  the  error  in 
(1],  it  is  beneficial  to  repeat  equation  (23)  of  [1],  l.e. 


CO  QO 


«»00 


5 


'  2  2 


where  Ax  -  -Xj  and  ^y'’y2”^l  *  (2.3)  I^x^^  ’^l)  point 

'‘2  2 

/  \  i-  i. 

j  on  the  large  scale  surface  having  slopes  and  is 

2  2  2 

illuminated  and  zero  if  it  is  shadowed.  The  function  f  ,  Is  defined  in  [11. 

PP 

~ 

is  the  mean  square  height  of  the  large  scale  surface,  and  i®  t;he 

normalized  autocorrelation  function  of  the  large  scale  surface  height.  It 
should  be  noted  that  the  left  hand  side  of  (2.3)  represents  the  Fourier  trans¬ 
form  of  the  <•>  term  (from  AxAy-  space  to  k  k  -space)  and  (2.3)  is 
supposed  to  show  where  the  transform  variables  appear  in  the  result.  Such  know¬ 
ledge  is  essential  to  accomplishing  the  convolution  of  (2.3)  with  the  small 
scale  surface  height  spectrum  because  this  convolution  determines  * 

That  is,  in  the  convolution  k  and  k  must  be  replaced  by  k  -k  and 

cx  oy  ox  X 

k^y-k^  where  k^  and  k^  are  the  new  variables  of  integration  (see  (32)  of 

[1]). 


If  R(0,(ti)  in  (2.3)  had  been  replaced  by  its  precise  definition,  as  given 

bv  (1),  the  problem  of  determining  where  k  and  k  appear  in  the  right 

cx  oy 

hand  side  of  (2.3)  would  have  been  correctly  solved.  Unfortunately,  this  was 
not  done.  Instead,  rt(G,(|>)  was  incorrectly  written  as  (1  +  C^)  ^  and,  through 
the  use  of  trigononetric  manipulations,  was  expressed  in  terms  of  k^^  and 

k^^  (see  equations  (24)  through  (28)  in  [1]).  This  development  failed  to  rec¬ 
ognize  that  does  not,  in  general,  depend  upon  the  transform  variables  in 

(2.3).  That  is,  if  one  changes  the  transform  variables  from  k  to  k  and 

ox  X 

from  k  to  k  ,  C  would  still  only  depend  on  k  =  -2k  sin  0  cos  4  and 
oy  y  o  ^  r  ox  o 

k^^  =  “2k^  sin  0  sin  (|>  .  The  primary  consequence  of  this  error  is  to  provide  an 
incorrect  formula  for  Rf*)  for  use  in  all  equations  following  (28)  in  [1]. 

As  will  be  shown,  this  error  fortunately  has  negligible  consequence  on  the 
numerical  results  presented  in  [1]. 
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2,3  Correct  Analysis 

The  error  Identified  and  explained  above  can  be  rectified  by  determining 
the  functional  dependence  of  the  right  hand  side  of  (2.1)  upon  the  large  scale 
slopes  and  since  they  are  replaced  by  the  transform  variables 

and  .  This  can  be  done  In  a  relatively  straightforward  manner  by  general¬ 

izing  Smith's  [3]  results  to  the  case  where  k^  is  at  an  angle  tt/2  -  4»  with 
respect  to  the  y-axis  rather  than  directly  along  the  y-axis.  Such  a  general¬ 
ization  leads  to  the  following; 


U(ctn0  -  Cj^^cosi,S  -  sin(t)) 


1  +  C 


(2. A) 


where  Is  given  by  (24) of  [1]  and  U(»)  is  the  unit  step  function  which 

is  one  If  the  argument  is  positive  and  zero  if  the  argument  is  negative.  Of 
particular  note  In  (2,4)  is  the  fact  that  and  appear  only  in  the 

argument  of  the  unit  step  function.  Substituting  (2.4)  into  (2.1)  yields  the 
correct  expression  for  R(6,<t))  ,  i.e. 


R(0,(f>) 


k  k 

U(ctn  0 - cos  (j) - ^  sin  (j)) 

o  B 

i  +  C„ 


(2.5) 


If  one  substitutes  k  =-2k  slnOcosd)  and  k  =  -  2k  sin  0  sin  d)  in  (2,5) 

such  as  required  in  the  determination  of  [0^°,]^  ,  then  R(0,4))  =  (1  +  C^)  ^ 

and  It  will  depend  only  on  the  angles  @  and  (j>  and  the  mean  square  slopes 
2  2 

and  of  the  large  scale  surface  (see  equations  (24)  and  (25)  of  [1]). 

However,  in  the  convolution  expression  for  unit  step  function 

must  be  retained.  That  is,  one  must  use  the  following  relationship; 


r  k  -k  k  -k  ' 

ox  X  oy  V 

k2k  cos  0  *  2k  cos  0 . 


Equation  (2.6)  leads  to  a  completely  different  Interpretation  of  the  ef¬ 
fects  of  large  scale  shadowing  on  ^*^pp*^j^  from  that  erroneously  presented  in 
[1],  According  to  (2.6),  shadowing  now  gives  rise  to  an  attenuation  factor 


(1+C  )  ^  which  is  common  to  both  [o  and  [a  .  This  result  is 

o  -  PP  0  -  PP  1 


o 

pp-  '0  -  ■  pp' 

merely  a  consequence  of  the  fact  that  small  areas  on  the  surface  capable  of 


producing  strong  Bragg  scatter,  i.e.  =  0  ,  are  as  equally  shadowed 

as  the  areas  properly  oriented  for  specular  reflection,  i.e.  “  ^ox^® 

=  ^oy^^  ’  statement  ;an  be  easily  verified  by  noting  that  for  both 

sets  of  the  above  values  of  and  >  U(ctn6  -  cos  <J)  -  sin  <t))  =  1 

in  equation  (2.4)  and  so  R(6,<|>)  ®  (1+C  )  ^  .  Since  both  the  I^nn'^O 

[0  °,1,  terms  suffer  the  same  attenuation  due  to  shadowing,  the  transition 
PP  1 

region  in  9  (where  f^pp«^Q  decreases  and  f'^pp*^!  becomes  predominant) 
is  independent  of  shadowing  effects.  Since  shadowing  results  from  the  slopes 
of  the  large  scale  surface  structure  and  since  the  large  scale  slopes  are  not 
the  important  surface  characteristic  in  determining  the  transition  region, 
this  result  demonstrates  that  the  theory  is  self  consistent. 

The  unit  step  function  in  (2.6)  serves  the  very  important  purpose  of  es¬ 
tablishing  the  limits  on  the  Integrals  in  the  convolutional  expression  for 

[0  °i]i  (see  (32)  of  [1]).  Hov’ever,  before  this  aspect  of  the  problem  is  con- 
PP  1 

sidered,  it  is  worthwhile  reviewing  the  physics  behind  the  reason  for  the  unit 
step  function  in  (2.6).  The  unit  step  function  appears  in  (2.6)  because  there 
is  a  certain  range  of  surface  slope  values  for  which  the  probability  of  having 
an  incident  ray  shadowed  is  identically  one  [3].  This  result  may  be  readily 
understood  by  referring  to  Figure  2,1.  In  (a),  the  slope  of  the  surface  is 
negative  at  the  point  of  Intersection  and  the  incident  ray  is  not  shadowed  in  a 
small  neighborhood  of  the  point.  In  (b),  the  surface  slope  is  equal  to  the 
slope  of  the  incident  ray  (ctnO)  and  the  ray  is  therefore  tangent  to  the  surface 
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. ...A.. . 


Figure  2-1 „  Diagrams  explaining  the  reason  for  the  unit  step 

function  in  the  shadowing  function.  Only  portions 
of  the  large  scale  surface  are  shown  since  the  small 
scale  structure  does  not  impact  the  shadowing. 


at  the  point  of  intetsection.  In  this  case  the  point  may  or  may  not  be  shadowed 
In  a  small  neighborhood  of  the  point,  depending  upon  the  surface  curvature.  In 
(c),  the  surface  slope  exceeds  the  slope  of  the  incident  ray  and  the  ray  is 
necessarily  shadowed  by  some  portion  of  the  surface  in  a  small  neighborhood 
of  the  point.  nus  all  points  on  the  surface  will  be  shadowed  if  their 
slope,  in  the  direction  of  the  incident  ray  projected  onto  mean  plane,  exceeds 
the  slope  of  the  incident  ray.  Stated  another  way,  the  probability  of  such 
an  event  is  one.  This  is  the  physical  reason  for  the  unit  step  function  in 
(2.4)  and  subsequent  equations  involving  the  shadowing  function. 

In  order  to  retain  the  physical  significance  of  the  «irilt  step  function, 
it  is  desirable  to  deal  with  a  particular  form  of  the  equation  for  (6,<j))  ]  , 

i.e.  equation  (40)  of  (1], 

oe»  a> 


T  (1  +  C  )  ^  i  J 

—00  —00 


cos  G  ^  +  k.  ,  2k  cos  6  C  +  k  ) 
^x  ox  ’  o  ^y  oy 


’  U(ctnO+^  cosd)  +  ^  sin  <j>)  F  ,  K  »  ~K  ) 
’x  ^  y  pp  x’  > 


•  exp 


.2  -I 


2C 


Zx 


2C 


Zy-* 


y.  y  k. 


(2.7) 


where  I  is  defined  in  [11.  Equation  (2.7)  expresses  the  convolutional 

broadening  of  the  spectrum  about  the  Bragg  wavenumbers  k^^  and  k^^  as  a 

direct  consequence  of  the  distribution  of  large  scale  slopes;  that  is, 

and  are  equivalent  to  and  *  It  should  be  pointed  out  that  the 

correct  expression  for  the  shadowing  function  has  been  used  in  (2.7).  The 

curve  in  the  C  C  -plane  separating  the  regions  where  the  step  function  is 
X  y 
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zero  and  one  is  a  straight  line,  l.e.  cos  (()  +  f;  sin  (()  ■  -  ctn6  .  Except 

X  y 

for  special  values  of  (|)  given  in  Table  1,  the  step  function  in  (2.7)  will 

consequently  give  rise  to  a  coupling  between  the  ^  and  Q  Integrals.  More 

X  y 

specifically,  the  lower  (upper)  limit  on  the  integral  in  (2.7)  is  given  by 


-  (ctnO  +  sin  <j)) 
cos  (|) 


(2.8) 


for  cos  0  >  0  (<0)  .  From  an  analytical  point  of  view  (2.8)  represents  an 
irritating  consequence  of  shadowing.  From  a  practical  standpoint,  the  restric¬ 
tions  Imposed  by  (2.8)  may  not  be  numerically  relevant  for  a  large  range  of 
incidence  angles  as  demonstrated  by  the  following  argument.  The  dominant  fac¬ 
tor  in  the  Integrand  in  (2.7)  is  the  slope  dependent  Gaussian  term  which  is 
equivalent  to  the  probability  density  function  for  the  large  scale  slopes.  For 
most  practical  purposes,  the  effect  of  this  term  is  to  truncate  the  range  of 
integration  in  (2.7)  to  about  a  ± 3-slgma  excursion  from  Cv  “  Cv  “  0  .  The 


±  3-slgma  excursions  for  C  and  ^  are 

X  y 


±3/^  and  ±3*/^. 


respec¬ 


tively,  and  substituting  these  values  in  the  unit  step  function  argument 
yields  the  following  requirement 


±  1  /  cos  (j)  ±  3  »  sin  (|)  >  -  ctn0 


(2.9) 


for  the  unit  step  function  to  be  unity.  If  0  <  <()  <  tt/2  ,  the  "worst  case" 

situation  occurs  when  C  =  -  zJ and  ^  =-3V^Cn^  or 

X  Jlx  y  Zy 


coset  +  3>/C^ 


3  /  cos  <|)  +  3  /  sin  (ji  <  ctn  0 


(2.10) 


The  worst  case  situation  occurs  when  the  left  hand  side  of  (2.9)  is  most  nega¬ 


tive.  Thus  for  numerical  purposes  and  all  values  of  9  satisfying  (2.10), 


one  can  replace  the  infinite  limits  in  (2.7)  by  and  ±3yCj^y 
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TABLE  I 


Integration  Limits  fur  Eqn.  (2.7)  and 
4*  a  Multiple  of  17/2 


4> 

intei 

llmi 

lower 

;  rat ion 

Lts 

upper 

5y  Intel 
lim: 

lower 

gratlon 

Lts 

upper 

0 

-ctn  G 

00 

«00 

00 

7r/2 

^  00 

00 

-  ctn  6 

oo 

TI 

—  00 

ctn  0 

—00 

00 

37T/2 

-  M  1 

«.  00 

00 

—00 

ctn  6 
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Conditions  similar  to  (2«10)  can  be  obtained  for  other  ranRes  of  (j)  In  a 
straightforward  manner. 

It  Is  Implicitly  assumed  that  the  upper  limit  on  0  resulting  from  (2.10) 

also  satisfies  the  optical  criterion  required  of  the  large  scale  surface*, 

2  2  2 

namely  that  4k^  cos  0  »  1  .  For  moderate  slopes  there  will,  generally  be 

a  gap  between  the  maximum  value  of  0  resulting  from  (2.10)  and  the  upper 

bound  on  0  resulting  from  the  large  scale  surface  optical  criterion.  In 

this  case,  one  must  necessarily  revert  to  the  more  exact  limits  such  as  given 

by  (2.8).  What  is  happening  in  this  situation  is  that  the  ±  3-sigma  support 

of  the  slope  density  function  in  (2.7)  is  overlapping  the  region  of  the 

5  5  -plane  where  the  unit  step  function  is  zeio.  In  fact,  if  one  could  go  to 
X  y 

the  grazing  incidence  limit  (0*tt/2),  it  is  readily  observed  from  Figure  2.2 

that  only  half  the  C  C  -plane  is  encompassed  by  the  integrals  in  (2.7).  Thus, 

X  y 

in  addition  to  the  attenuation  of  near  grazing  incidence  due  to  the 

(1+C^)^  factor,  there  is  another  reduction  factor  resulting  from  the  shadowing 
of  points  on  the  large  scale  surface  having  positive  slopes  (see  Figure  2.1). 

It  should  be  remembered  that  although  the  impact  of  the  unit  step  function  upon 
the  limits  of  the  integrals  in  (2„7)  is  somewhat  involved,  it  is  a  direct  con¬ 
sequence  of  the  rather  simple  slope-shadowing  limitations  explained  in  Figure 
2.1. 

The  analyslc  pi  .'sented  above  does  not  alter  the  general  composite  surface 
scattering  theory  set  forth  in  [1],  It  does,  however,  correct  and  expand  the 
theory  in  [1]  a/i  it  relates  to  the  effects  of  shadowing  upon  large  angle  of 
incldance  backf/c atcering  frou  a  randomly  rough  surface.  This  additional  analy¬ 
sis  was  necessitated  by  the  use  of  an  incorrect  functional  form  for  the  shadowing 
fuTiCtion  in  fl). 

*Thi3  criterion  has  also  been  called  the-  stationary  phase  approximation  [4] 
and  the  dfieu  phase  modulation  condition  [51. 


2, A  Discussion  of  Numerical  Example 


The  use  of  an  incorrect  shadowing  function  in  Sections  I  through  IV  oi 
[1]  was  primarily  a  sin  of  omission.  That  is,  since  nearly  all  the  results 
in  Section  I  through  IV  c f  [1]  were  formal  In  nature,  it  was  not  necessary 
to  use  the  erroneous  shadowing  function  given  by  (28)  of  [1].  However,  the 
numerical  example  comprising  Section  V  of  [1]  does  require  attention  in  order 
to  properly  account  for  the  effects  of  large  scale  shadowing.  In  the  example 
presented  in  [1],  R  was  set  equal  to  unity  because  the  large  scale  slopes 
were  relatively  small.  According  to  the  analysis  presented  above,  this  step 
is  not  justified  for  all  angles  of  incidence,  Tlie  correct  effect  of  shadowing 
is  addressed  below. 

2  2 

As  in  Section  V  of  [1],  it  is  assumed  that  w  (0,C)  and  the 


surface  height  spectrum  is  Isotropic,  l.e. 


.2  ^2  ^  2  . 

_ _  *«.x  “  ^’«,y  "  ^«,t 

S(kx.ky)  -  s(t/k^  ky  )  .  Substituting  k  cos  a  »  ,  k  sin  a  ■  and 

dk  dk  ■  kdkda  in  (32)  of  [1]  and  using  (2.6)  above  for  the  shadowing  func- 
X  y 

tion  yields; 


'Vp-’i 


u  cos^9(l  +  C^> 


2T[  00 

r  f  S(k)u(ctnO  +  tan  0  + 

J  •/  o 

o  k , 
d 


•  exp  [AkJ^  sln^  0  +  +  Akk^  sinO  cos(a-(|))  ]  j Ak^  cos^  0  |  kdkda  (2.11) 

The  unit  step  function  is  unity  whenever  kcos(a-(|))  >  -2k^/8in0  .  Since  the 

minimum  value  of  co8(oi— ({))  is  -1,  this  inequality  will  be  satisfied  for  all 

a  if  k  is  restricted  to  less  than  2k  /sin0  .  Since  the  range  of  k  from 

o 

2k  /sind  to  “  for  cos(a-(())  >  0  contributes  little  to  the  integral  in  (2.11) 
o  _____ 

2 

(because  is  smell),  the  k-integral  limits  can  be  approximated  by 
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[kj,2k^/ain  0]  .  In  this  range  of  k  ,  the  step  function  Is  unity  and  the  a- 
integratlon  can  be  accomplished  with  the  following  result; 


4k^  r  ^,(0,0) 

[o  °,]  :: - 2__££'' 

pp’  I  “T  2 

cos  0(1+C^ 


2k  /sine 

o 


.0)  r 
+c )  y 


s(k)  I  I  -AiiSLL- 

'  o{  — y 


exp  < 


^  (k--2k^  8ln0)^ 

27~2  2  . 

'"‘o  ® 


k  sin  0 


k  dk 


(2.12) 


where  I  (*)  Is  a  Bessel  function  of  the  second  kind.  The  right  hand  side  of 
o 

(2,12)  may  now  be  compared  with  the  [a  part  in  (44)  of  [1],  i.e. 

PP  * 


^0  r 

7^  2  o  J 

hr.  ®  k. 


S(k)  I 


k  sin  6 


ko  cos^  0 


exp  <  - 


(k-2k  sine) 

0 _ 

2'ri  2q 
4ko  cos  0 


k  sin  0 


k  dk 


(2.13) 


The  obvious  differences  are  the  factor  (1+C^)~^  and  the  finite  upper  limit 

on  the  integral  in  (2.12).  Figure  2.3  illustrates  how  (1+C^)  ^  varies  with 
2 

0  for  =  0.0224  which  was  the  value  of  mean  square  slope  used  to  construct 

Figures  3  and  4  of  [1],  Of  particular  note  in  the  plot  of  (1+C^)  ^  is  the 
fact  it  does  not  start  to  decrease  until  0  exceeds  85“ ;  at  87.5" 

(1  +  C  )  ^  =  -2.5  dB.  The  value  of  0  =  87,5"  is  the  point  at  which  4k^Cj  cos^0~lO 
and,  consequently,  represents  the  approximate  limit  of  the  large  scale  theory. 

That  is,  for  0  >  87.5"  the  analysis  of  the  scattering  from  the  large  scale 
structure  on  the  surface  can  no  longer  be  accomplished  using  optical  techniques. 
The  effect  of  the  finite  upper  limit  on  the  Integral  in  (2,12)  is  much  less 
dramatic.  In  fact  it  the  4kJ^J^j^cos  0  ^  10  criterion  is  Ignored  and  the  limit 
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of  0  ■  90“  is  taken  in  (L.12),  it  can  be  sbtvm'^  that  the  finite  upper  limit 
gives  rise  to  a  3  dB  reduction  from  the  value  o^  (2.13)  at  6  ■  90“ .  This  oame 
conclusion  can  also  be  obtained  from  (2.7)  and  it  is  a  direct  consequence  of 
the  fact  that  all  positive  slopes  are  excluded  from  considoration  at  0  ■  90“ 

(see  Figure  2.1).  Within  the  range  of  validity  of  the  optical  criterion,  the 
finite  upper  limit  on  the  Integral  in  (2.12)  gives  rise  to  much  less  attenua¬ 
tion  than  the  (1+C  )  factoro 

o 

In  view  of  the  above  analysis,  it  is  concluded  that  the  numerical  results 
presented  in  Section  V  of  [1]  are  correct  as  they  were  presented.  Determina¬ 
tion  of  the  spectral  division  wavenumber  is  in  no  way  altered  by  the  inclu¬ 
sion  of  the  correct  shadowing  function.  The  curves  shown  in  Figures  3  through 
6  of  [1]  are  correct  because  they  do  not  encompass  the  range  of  0  where  shadow¬ 
ing  is  important  (0~85“).  As  noted  above,  the  onset  of  shadowing  effects 
occurs  approximately  where  the  optical  criterion  (4k  ^  cos^  0  -  10)  is  vlo- 

O  X* 

lated,  i.Bj  6  ~  87.5“  ,  for  the  numerical  example  presented  in  [1],  However, 
for  larger  slopes  shadowing  must  be  considered  since  it  will  cause  a  signifi¬ 
cant  reduction  in  '^pp«  near  grazing  incidence. 

2 . 5  Summary 

In  the  analysis  presented  in  [1],  it  was  correctly  demonstrated  how  one 
includes  the  effects  of  large  scale  surface  feature  shadowing  on  the  backscat- 
tering  cross  section  of  a  composite  surface  for  large  angles  of  incidence. 
Unfortunately,  an  incorrect  form  for  the  shadowing  function  was  used  in  [1] 

which  led  to  the  erroneous  evaluation  of  the  impact  of  shadowing  upon  large 
angle  of  incidence  scattering.  In  this  section,  the  correct  form  of  the  shadowing 

*To  sliow  this  one  can  apply  Laplace’s  method  to  asymptotically  evaluate  (2.12) 
as  cosO->-0.  However,  it  must  be  remembered  that  the  maximum  of  the  integrand 
occurs  at  the  upper  limit  of  the  integrand  as  cos0->O  and  this  Impacts  the 
evaluation  of  the  Integra]  16]. 
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function  has  been  presented  and  Included  in  the  formulas  obtained  In  [1]. 
Particular  emphasis  has  been  placed  upon  the  physical  significance  of  shadow** 
Ing  as  it  effects  large  angle  of  incidence  scattering.  It  has  been  shown  that 
shadowing  leads  to  multiplication  of  both 

tor  (1+C  which  causes  (o  and,  thus,  o  to  go  to  zero  near 

o  pp'  1  pp' 

grazing  Incidence  for  a  perfectly  conducting,  randomly  rough,  composite  sur¬ 
face.  Furthermore,  there  is  another  effect  which  leads  to  an  additional  3  dB 
attenuation  at  grazing  incidence  (0^-90“)*.  This  effect  results  from  those 
slopes  which,  with  probability  one,  will  cause  the  point  on  the  surface  having 
these  slopes  to  be  shadowed.  At  grazing  incidence,  all  positive  slopes  are  in 
this  class. 

A  reevaluation  of  the  numerical  results  presented  in  [1]  revealed  that 
use  of  the  coi , ect  shadowing  function  did  not  alter  any  of  the  results  rela¬ 
tive  to  the  choice  of  the  spectral  dividing  wavenumber  k^.  Furthermore,  none 
of  the  curves  presented  in  [1]  were  affected  because  they  only  encompass  the 
range  of  0  from  0  to  70“  and  the  effects  of  shadowing  were  present  for  6-85“. 
Also,  techniques  were  presented  relative  to  overcoming  some  of  the  analytical 
difficulties  resulting  from  the  use  of  the  correct  shadowing  function. 


*It  is  reemphasized  that  exact  grazing  incidence  cannot  be  addressed  by  this 
theory  because  the  optical  criterion  assumed  of  the  large  scale  surface  fea¬ 
tures  is  violated.  The  -  3  dB  figure  is  significant  only  in  its  magnitude  rela¬ 
tive  to  the  effect  of  the  (1+C  factor. 

o 
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3.0  SHADOWING  BY  NON-GAUSSIAN  RANDOM  SURFACES 


3.1  Background 

Shadowing  of  random  surfaces  was  originally  introduced  [1]  as  an  ad  hoc 
correction  to  the  results  provided  by  physical  or  geometrical  optics  approxi¬ 
mate  theories  of  rough  surface  scattering.  Sancer  [2]  subsequently  demonstr.i— 
ted  how  shadowing  could  be  rigorously  accounted  for  in  the  optical  limit  for 
random  surfaces.  Furthermore,  he  showed  that  previously  derived  expressions 
for  the  effects  of  shadowing  based  upon  purely  geometrical  considerations  [3, 

4]  were  directly  applicable.  Using  Sancer's  results.  Brown  [5,6]  showed  how 
shadowing  could  be  rigorously  included  in  a  formulation  for  scattering  from 
random  surfaces  characterized  by  many  scales  of  roughness,  i.e.  composite 
rough  surfaces. 

While  shadowing  theory  is  reasonably  mature.  It  has  only  been  applied  to 
jointly  Gaussian  random  surfaces.  The  Gaussian  results  are  probably  adequate 
for  the  ocean  but  they  are  questionable  for  terrain  and  completely  inadequate 
for  sea  ice  fields.  For  sea  ice,  water  first  fills  all  surface  depressions 
below  mean  sea  level  and  then  freezes.  This  eliminates  all  surface  height 
excursions  below  mean  sea  level  and  the  probability  density  function  of  the 
surface  roughness  is  clearly  non-Gausslan.  For  these  reasons,  it  is  important 
to  extend  shadowing  theory  to  the  point  where  it  can  easily  accommodate  non- 
Gausslan  surface  statistics;  such  is  the  purpose  of  this  section. 

3.2  Analysis 

The  special  case  of  backscattering  is  chosen  to  illustrate  the  approach; 
this  minimizes  some  of  the  conceptual  details  associated  with  the  more  general 
bistatic  case.  It  turns  out  that  the  extension  of  the  results  to  the  blstatlc 
geometry  can  be  accomplished  almost  by  inspection.  The  analysis  presented  by 
Smith  [4]  is  general  to  a  point  in  the  development;  however,  there  are  a  number 
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of  Integrations  which  must  be  accompllslied  In  order  to  arrive  at  the  final 
expression  for  the  shadowing  function.  In  the  case  of  a  jointly  Gaussian 
surface  these  integrals  can  be  performed  and  a  closed-form  result  is  obtained 
for  the  shadowing  function.  For  non-Gausslan  surfaces,  the  required  integra¬ 
tions  appear  to  be,  at  best,  formidable.  The  purpose  of  this  section  is  to 

show  that  if  the  height  and  slopes  of  the  surface  are  independent  random  vari¬ 

ables  then  the  final  expression  for  the  shadowing  function  is  drastically  sim^ 
plifled. 

For  the  reader's  convenience  Smith's  notation  will  be  employed  in  this 
section  and  his  Figure  1  Is  essentially  repeated  here  as  our  Figure  3.1.  There 
are  three  critical  relationships  from  Smith's  paper  [4]  which  are  required. 

If  S(F,0)  is  the  probability  that  no  part  of  the  surface  will  intersect  the 
incident  ray  (at  an  angle  0  with  respect  to  the  normal  to  the  mean  flat  sur¬ 
face)  on  its  way  to  point  F  on  the  surface  then  S(F,0)  is  given  by 


S(F,0)  =  h(y  -q^)  exp 


CO 

1^  -  J*  g(T)dT 


(3.1) 


where  h(»)  is  the  unit  step  function,  p  =  ctn  0  ,  q^  is  the  slope  of  the 
surface  in  the  y-direction  at  F  ,  and  g(T)AT  is  the  conditional  probability 
that  the  surface  will  intersect  the  incident  ray  in  the  interval  (t,T+At) 
given  that  it  does  not  intersect  the  ray  in  (0,T)  .  The  function  g(T)AT  Is 
determined  by  the  behavior  of  P2(C,q|F,T)  which  is  the  joint  probability  of 
the  height  C  and  y-slope  at  the  point  (x  =  0,y  =  T)  conditioned  upon  the 
height  (Cq)  and  y-slope  (q^)  at  point  F  ;  in  particular. 
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Figure  3-1.  Shadowing  geometry.  The  incident  ray  lies  in 
the  X  =  0  plane  and  the  slopes  of  the  surface 
the  point  F  are  35/3x  =  p  and  3i^/3y  =  q 


o 

f 


(q  -y)P3(f,.qlF,T) 


dq 


g(T)AT 


®  Cq+WT 

J  dq  J  P3a.qlF, 


•At 


T)dC 


(3.2) 


The  average  oE  S(F,0)  over  all  surface  heights  with  ■  0  and  q^  ■  -  1/y 
is  the  desired  shadowing  function  R(0)  ,  i.e.  the  probability  that  a  back- 

scattering  specular  point  on  the  surface  will  not  be  shadowed. 

Smith  proceeded  to  evaluate  (3.2)  in  the  Gaussian  case  by  assuming  that 
the  heights  and  slopes  at  F  were  uncorrelated  with  those  at  y  =  T  .  Here 
it  will  be  assumed  that  decorrelation  implies  statistically  independence 
and  that  tlie  height  and  slopes  are  independent;  thus, 

P3(C,qlF.T)  P^(^)P2(q)  (3.3) 


where  Pj^(?)  is  the  height  probability  density  function  and 


CO 

■I 


(3.4) 


where  P22(p.q)  Is  the  joint  probability  density  function  of  the  x  and  y 
slopes.  Substituting  (3.3)  into  (3.2)  yields  the  following 


g(T)  = 


r  ?,(£  +yT) 
1  o 


I 


Co+UT 


(3.5) 


Pl(C)dC 


where 


•This  is  not  true  in  general  so  the  following  analysis  applies  to  a  restricted 
class  of  non-Gaussian  surfaces. 
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r  -  I  (q  -M)P,(q)dq 


(3.6) 


The  denominator  of  (3.5)  is  recognized  as  the  distribution  function  for  C 
evaluated  at  +1^^  •  i.e.  Fj^(^^+Ijt)  .  Also,  Pj^(C^+l^ii)  in  (3.5)  is 

equal  to  the  derivative  of  the  distribution  function  evaluated  at  C  =  r.  +  MT 

O 

Consequently,  (3.5)  becomes  g(T)  -  tr/Pj^(CQ +yT)]dFj^(i;^ +lJT)/d(^^ +VIT)  . 

Substituting  this  result  in  (3.1),  making  the  change  of  variable  +yT 

o 

and  noting  that  dF^^/Fj^  -  d(lln  Fj^)  yields 


OO 

S(0,F)  *  h(y-q^)  exp  ^^-T/y  J  d[J,nFj^(n)  ] ^ 


(3.7) 


where  An  denotes  the  natural  logarithm.  Since  Fj^(")  =  1  ,  (3.7)  reduces  to 


the  following; 


S(e,F)  =  h(p-q^)  F^  (C^)' 


(3.8) 


Remembering  that  ^2^(^q)  *  average  of  (3.8)  over  all  values 

of  simplifies  to 


00 

S(p„,q„.6)  -  h(M-,„)  J  dFj(e„) 


(3.9) 


h(y  -  q  ) 
s(Po.qo»9)  =  r/y  +  1 


(3.10) 


since  Fj^(-“>)  =  0  and  F/y  +  1  ^  0  .  With  qQ  “  ”  l/y  *  the  shadowing  func¬ 

tion  appropriate  for  backscatter  reduces  to  the  following  simple  expression; 


V 


00 


R(0) 


(q/U 


i)i’2(‘i)dq 


') 


-1 


(3.11) 


where,  in  summary,  y  “  ctn  6  ,  6  Is  the  Incidence  angle,  and  P2  (q)  is 
the  probability  density  function  of  the  slopes  of  the  surface  in  the  plane  of 
Incidence  defined  by  the  Incident  ray  and  the  normal  to  the  mean  surface.  It 
is  interesting  to  note  from  (3.11)  that  for  normal  incidence  (0“O)  R(0)  ■  1 

whereas  at  grazing  incidence  (0  =  tv/2)  R(lf/2)  ■  0  since  li  “  0  and 

00 

J  qP2(q)dq  >  0  .  Thus,  these  basic  properties  of  the  shadowing  function  are 
o 

independent  of  the  detailed  properties  of  the  slope  density  function.  Cfne  can 
easily  verify  that  (3.11)  is  identical  to  the  results  obtained  by  Smith  for 

the  special  case  of  a  jointly  Gaussian  surface. 

The  form  of  (3.10)  compared  to  Smith's  results,  i.e.  (23)  of  [4],  suggests 
that  the  above  result  can  be  directly  translated  to  the  blstatlc  case  and,  in¬ 
deed,  this  is  the  case.  For  the  blstatlc  case,  a  generalization  of  Sancer's 
[2]  results  will  be  given.  It  should  be  noted  that  the  inequalities  involv¬ 
ing  the  angles  of  incidence  (®q)  scattering  (0)  ,  just  prior  to  Sancer's 

equations  (49),  (50),  (54)  and  (55)  should  be  reversed.  With  y  * ctn  0  and 
y^  = ctn  0^  ,  Sancer's  results  are  easily  generalized  by  replacing  his  by 

r{y^)/y^  and  C2  by  r(y)/y  . 


3. 3  Example 

To  Illustrate  the  above  results,  the  backscattering  shadowing  function 
R(0)  will  be  determined  for  the  exponential  joint  slope  density  function  in¬ 
troduced  by  Barrick  [7],  i..e. 

P22(p»q)  [“  V^0(p^  +  q^)/w^  ]  (3.12) 

TTW  ^ 
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where  w  l8  the  mean  square  slope  of  the  surface  roughness.  P22(p»q)  in 
'3.12)  represents  a  surface  whose  roughness  is  isotropic,  w  /2  “  ‘^p  >  “  ^q  ^  > 
with  statistically  dependent  slopes.  That  is,  the  joint  density  function  can¬ 
not  be  expressed  as  a  product  of  the  marginal  or  individual  densities.  The 
calculation  of  the  marginal  density  »  ““log  (3. A),  is  reasonabiy  straight¬ 

forward  and  the  result  is  as  follows; 

P2(q)  -  hi  Ki  (/6  |q|/w)  (3.13) 

TTW 

where  Kj^(*)  is  one  of  the  modified  Bessel  functions  of  order  one  t8].  It 
■*8  interesting  to  compare  this  density  with  a  Gaussian,  i.e. 

®  ~ —  exp(-q^/2w^) 

and  this  is  done  in  Figure  3.2  where  the  normalized  densities  P2(q)w  are 
plotted  as  a  function  of  the  normalized  slope  q/w  .  It  should  be  noted  from 
f*  'lots  in  Figure  3.2  that  the  "exponential"  density  shows  a  much  greater 
prc  'bility  of  occurrence  of  small  slopes  than  the  Gaussian.  This  result  is 
in  agreement  with  one  intuitive  approach  for  generating  a  surface  character¬ 
ized  bj  (3.12),  e.g.  one  strongly  filters  all  surface  height  excursions  below 
a  ce'  in  level  to  eliminate  the  possibility  of  large  negative  height  excur¬ 
sions.  This  process  Increases  the  probability  of  small  slopes  at  the  expense 
of  the  large  slopes. 

Substituting  (3.13)  in  (3.11)  and  using  tabulated  integrals  of  Bessel 
functions  given  in  [9],  the  following  closed-form  result  is  obtained  for  the 
backscatterlng  shadow  function  R(9)  ; 

R(e)  »  K2(x)  +  I  +  f  [^K^(x)L^(x)  +  L^(x)K^(x)]  ^  (3.14) 
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where  x  -  (>/6/w)ctn  0  and  the  ,  n=0  or  1  ,  ayribol  denote.,  the  modified 

Struve  functions  [8].  Using  asymptotic  forms  for  the  special  functions  in 
(3.1A),  it  may  be  readily  verified  that  R(0)  -  l(x^)  and  R(it/2)  >«  0(x«0)  . 
The  modified  Struve  functions  may  be  computed  from  tables  given  in  [8]  for 
X  -  5  and  by  a  power  series  for  smaller  arguments. 

Figure  3.3  compares  (3.1A)  and  the  shadowing  function  for  a  Gaussian  func 
tion  obtained  by  Smith  [A]  for  a  range  of  rms  slopes.  The  shadowing  function 
for  the  exponential  Joint  slope  density  is  larger  because  the  marginal  slope 
density  given  by  (3.13)  exhibits  less  likelihood  for  large  slopes  than  the  cor 
responding  Gaussian  density.  That  is,  the  larger  slopes  are  the  source  of 
more  significant  shadowing. 

3. A  Summary 

The  shadowing  theory  developed  by  Smith  [A],  while  sufficiently  general 
to  deal  with  any  joint  slope  density  function,  involves  what  appears  to  be  a 
number  of  rather  complicated  Integrals.  Under  the  assumption  that  the  sur¬ 
face  height  is  statistically  independent  of  the  surface  slopes,  it  is  shown 
that  Smith's  theory  can  be  reduced  to  a  single  Integration  involving  the  mar¬ 
ginal  density  function  for  the  slopes  in  the  plane  of  incidence.  Using  this 
result  but  without  regard  to  the  specific  form  of  the  marginal  density  func¬ 
tion,  it  can  be  shown  that  the  backscattering  shadowing  function  is  unity  at 
normal  Incidence  and  zero  at  grazing  incidence.  Because  the  final  result  in¬ 
volves  an  integration  or  smoothing  process,  it  is  amenable  to  the  use  of 
histogram  data  for  the  marginal  slope  density. 

This  theory  is  applied  to  an  exponential  joint  slope  density  represent¬ 
ing  an  isotropic  surface  for  which  the  slopes  are  not  statistically  indepen¬ 
dent.  The  backscattering  shadowing  function  for  the  exponential  and  Gaussian 
joint  slope  densities  are  compared  and  it  is  found  that  the  Gaussian  surface 
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Figure  3-3.  A  comparison  of  the  backseat coring  shadowing  function  for  the  exponential 
and  Gaussian  joint  slope  densities  for  a  range  of  rms  slopes  (w). 


produces  stronger  shadowing.  This  result  is  found  to  be  a  consequence  of  the 

greater  likelihood  of  large  slopes  with  the  (laussian  density. 
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4.0  BISTATIC  SCATTERINC  FROM  LOSSY  RANDOM  SURFACES 


4.1  BnckRround 

Prior  to  the  iMld-1960'8,  electronagnetlc  scattering  from  randomly  rough 
surfaces  was  modeled  using  either  perturbation  theory  or  physical  optics  [1], 
First  order  perturbation  theory  appeared  to  do  a  reasonable  job  of  analytically 
describing  the  scattering  process  when  the  surface  roughness  was  small  in  terms 
of  the  electromagnetic  wavelength  and  multiple  scattering  was  negligible. 
Physical  optics  produced  meaningful  results  in  and  about  the  specular  scatter¬ 
ing  direction  wlien  tlie  surface  exIUbited  very  large  but  smoothly  undulating 
height  variations.  Unfortunately,  there  were  numerous  attempts  to  apply  these 
theories  to  situations  where  the  Implicit  assumptions  in  the  models  were  vio¬ 
lated.  These  attempts  usually  assumed  some  surface  parameter  such  that  the 
scattering  measurements  and  the  "model”  were  brought  into  agreement.  However, 
it  was  very  quickly  recognized  that  these  attempts  were  highly  suspect  because 
of  their  failure  to  satisfy  certain  fundamental  principles. 

As  more  and  more  rough  surface  microwave  scattering  measurements  were 
acquired,  It  became  obvious  that  neither  first  order  perturbation  theory  nor 
physical  optics  were  individually  adequate  for  all  angles  of  incidence  and 
scattering.  Conversely,  it  appeared  that  physical  optics  seemed  to  do  a  good 
modeling  job  near  the  specular  scattering  direction  while  first  order  pertur¬ 
bation  theory  was  reasonably  accurate  for  all  other  scattering  angles.  Almost 
simultaneously,  researchers  in  the  U.S.  |2]  and  the  U.  S,  S.R,  [  3,4 ]  begf.n  to 
advocate  the  combining  of  these  two  diverse  theories  in  what  was  later  to  be 
called  the  composite  surface  scattering  model.  In  this  model,  the  surface 
was  considered  to  be  made  up  of  both  large  and  small  scale  surface  features 
(height  and  spatial  wavelength)  relative  to  the  electromagnetic  wavelength, 

.  The  large  scale  surface  features  were  considered  to  be  responsible  for 
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the  physical  optics-like  scattcrlni>  near  the  specular  direction.  The  small 
scale  surface  structure  gave  rise  to  a  perturhatlon  field  (to  the  optical 
field)  which  was  the  dominant  scattering  mechanism  away  from  the  specular 
scattering  direction.  The  interaction  between  the  optical  and  first  order  per¬ 
turbation  fields  was  assumed  to  be  totally  dependent  upon  the  tilting  of  the 
small  surface  structure  by  the  larger  gently  undulating  features  [5], 

More  recently,  rigorous  first  order  boundary  perturbation  theory  has  been 
applied  to  the  problem  of  backscattering  from  a  perfectly  conducting,  Gaussian 
distributed  rough  surface  [6],  The  results  of  this  analysis  Indicated  that 
much  of  the  original  work  on  this  problem  could  be  rigorously  justified. 
Furthermore,  additional  insight  was  gained  in  regard  to  such  asp-'cts  of  the 
problem  as  shadowing  (see  Section  3  of  this  report),  spectral  dichoto’^y,  and 
the  tilting  interpretation.  A  logical  extension  of  this  latter  theory  encom¬ 
passes  blstatlc  scattering  from  a  lossy  dielectric  surface.  The  purpose  of 
this  section  is  to  present  the  details  associated  with  such  an  extension. 

Before  the  details  of  scattering  from  a  composite  dielectric  surface  are 
presented,  it  is  illuminating  to  consider  two  mucli  simpler  cases.  The  first 
is  backscattering  from  a  dielectric  surface  having  only  a  small  scale  rough¬ 
ness  while  the  second  case  addresses  blstatic  scattering.  The  advantages  of 
this  approach  are  that  it  leads  to  familiarity  with  the  perturbation  technique 
and  it  sets  forth  the  principles  that  will  be  used  for  the  composite  surface. 

4.2  Backscattering  From  A  Dielectric  Surface  With  Small  Scale  Roughness 

The  geometry  for  this  problem  is  shown  in  Figure  4-1,  Tlie  mean  or  aver¬ 
age  surface  is  the  z  =  0  plane;  the  random  roughness  superposed  upon 

s 

this  plane  is  positive  for  ?  >0  and  negative  for  <  0  .  Below  the  rough 

s  s 

surface  (z <  C  ),  the  relative  dielectric  constant  of  the  medium  is  e  and 
s  r 

the  relative  magnetic  permeability  is  taken  to  be  the  same  as  for  free  space. 
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for  bnckacattering  from  a  randomly  rough 
avlnc  only  small  scale  roughness  C 


i.e.  ^  •  Above  the  rough  surface  (z  >  ,  the  mcdltini  Is  free  space, 

i.e.  e  “1  and  y  =  1  . 
r  '^r 

Provided  that  the  roughness  is  small  with  respect  to  the  electromagnetic 

2  2  2 

wavelength  ,  e.g.  4k^  Cg  1  wlicre  =  2tt/X^  and  Cg  is  the  mean 
square  height  of  the  roughness,  the  scattered  field  can  be  expressed  as 

follows; 

E  =!  6  E  +  6  E  (4.1) 

where  6  E  is  the  field  scattered  by  a  surface  having  no  roughness  (the  zeroth 
order  perturbation)  and  6*E  is  the  scattered  field  wliich  depends  on  the  rough¬ 
ness  to  first  order  only  (the  first  order  perturbation).  The  primary  assump- 

2  3 

tion  in  (4.1)  is  that  higher  order  terms  such  as  0(C  )  ,  0(C  )  ,  etc. ,  are 

s  s 

negligible.  The  zeroth  order  perturbation  field  <S”e  is  trivially  determined 
since  it  is  Just  the  field  reflected  by  an  Infinite,  flat  dielectric  interface. 
Both  Mltzner  [7]  and  Burrows  [8]  have  obtained  particularly  useful  expressions 
for  6  E  .  The  Mitzner  result  is  more  straightforward  but  it  is  restricted  to 
small  roughness  perturbations  superposed  on  a  flat  plane.  Burrows*  solution 
for  o*E  is  somewhat  more  complicated  but  it  is  more  general  in  that  the  unper¬ 
turbed  surface  need  not  be  planar  or  even  deterministic.  For  small  scale  rough¬ 
ness  on  a  plane,  the  Burrows  formulation  requires  a  bit  more  effort  in  computing 
6^E  than  Mltzner* s  result.  For  a  composite  surface,  only  the  Burrows  result 
is  sufficiently  general  to  address  this  problem. 

At  first  glance,  the  Burrows  expression  for  o  E  appears  to  be  somewhat 
cumbersome  and  confusing.  However,  if  it  is  realized  that  the  result  is  obtained 


®Mltzner*s  result  can  actually  be  applied  to  any  unperturbed  surface  for  which 
the  wave  equation  is  separable.  For  the  problem  considered  here,  the  surface 
is  restricted  to  a  plane. 
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from  an  application  of  reciprocity  then  the  notation  becomes  more  meaningful. 
Basically,  one  deals  with  two  incident  electric  fields  of  the  form 

e  *^i  ~  ®xp(  -  j  k^»r)  (4.2a) 

and 

=  E^'  e*  E^'  =  exp(  -  j  •  r)  (4.2b) 

where  the  primed  field  may  have  a  different  polarization  and  direction  of  inci¬ 
dence  than  the  unprimed  field.  Burrows'  expression  for  the  first  order  per- 
turbation  (electric)  field  scattered  in  the  direction  -  and  polarized  in 
the  e'  direction  Is  as  follows  [8]; 


6‘E.e'  = 


o  o  J 


where  it  is  assumed  that  6*E»e'  is  measured  in  the  far-field  of  the  rough 
surface.  The  distance  R  is  measured  from  the  origin  of  the  reference  co¬ 
ordinate  system  on  the  mean  surface  to  the  point  of  observation  or  measurement 
of  6*E*e'  and  is  the  permittivity  of  free  space.  The  fields  E'  ,  H'  , 

D'  and  B'  are  the  fields  on  the  unperturbed  surface  (S^)  due  to  the  primed 
Incident  field  while  AE  ,  AH  ,  Al)  and  AB  are  the  discontinuities  in  the 
fields  on  the  unperturbed  surface  to  the  unprimed  incident  field. 

Thus,  to  determine  tlie  e  -polarized  component  of  the  scattered  first  order 

A 

perturbation  field  in  the  general  direction  one  merely  sets  e'  = e^  and 

k^  =  -  k^  in  the  expression  for  the  primed  incident  field  E^  ,  computes  the 

->■  '> 

resulting  fields  E'  ,  H'  ,  D'  ,  B'  on  the  unperturbed  surface  S  ,  and  sub- 

o 

stltutes  these  results  in  (4.3). 

A  A  A  A 

For  backscatterlng  k^  =  -k^  so  according  to  the  above  receipe,  k^'  =•  k^ 
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which  In  the  coordinates  of  the  geometry  shown  in  Figure  A-1  is  as  follows; 

“  k^(  -  sin  0  cos  <|>  X  -  sin  0  sin  4)  y  “  cosGz)  (4.4) 

Contrary  to  previous  analyses  [1],  the  direction  of  incidence  specified  by 
the  angle  4>  should  not,  at  this  point  in  the  development,  be  arbitrarily  set 
to  some  convenient  value  such  as  0  or  7t/2  .  The  reason  for  this  is  that  the 
surface  may  have  anisotropic  roughness  and  the  orientation  of  the  x  and  y 
axes  of  the  reference  coordinate  system  should  be  fixed  relative  to  this  sur¬ 
face  characteristic  and  not  the  direction  of  incidence. 

Since  the  fields  inside  the  surface  integral  in  (4.3)  are  the  fields 

induced  on  the  infinite  planar  dielectric  surface  S  ,  it  is  convenient  to 

o 

further  categorize  the  problem  according  to  "’•■e  polarization  of  the  incident 
fields.  For  both  and  horizontally  polarized,  e  and  e*  are 

orthogonal  to  the  plane  formed  by  the  unit  vectors  and  *  z  where  n 

is  the  normal  to  the  mean  or  unperturbed  surface.  In  this  case  both  e  and 
e’  are  totally  tangential  to  the  mean  plane.  For  both  E^  and  E^'  vcrti- 
cally  polarized,  e  and  e'  are  parallel  to  tl*v.  plane  formed  by  and 

A  A 

n  =  z  . 

4.2.1  Horizontal  Polarization 

A 

When  e  and  e’  are  tangential  to  the  mean  or  unperturbed  surface,  the 
fields  E  ,  E’  ,  D  and  D'  (when  evaluated  on  S^)  are  entirely  tangential  to 
.  Since  the  tangential  component  of  the  electric  field  is  continuous  across 
an  Interface,  AE  =  0  .  Furthermore,  since  there  is  no  change  in  across 
the  boundary  and  the  lower  medium  is  assumed  not  to  be  perfectly  conducting, 

AB  and  AH  are  both  zero  on  the  Interface.  Consequently,  (4.3)  reduces  to 
the  following; 
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. : : . 

-’^1  .  IHtWHWWIWVF  WHii>iWtv«v.j-  »•  .  i»^.iii,,.— ..  w.  .  . 

-  k  ^  exp(  -  jk  R)  r  ^ 

*  -  °4.RE  e  '  1 

o  o  J 

S 

o 

i 

.1 

j 

(4.5)  1 

1 

The  total  E* -field  on  the  surface  due  to  E^^'  is  given  by 

1 

•  i 

(1+Rj^)  E^  (S^)  =  E^(1+Rj^)  exp 

i 

(4.6)  1 

0 

i 

1 

a 

a 

j 

where  R^  is  the  Fresnel  (field)  reflection  coefficient  for  horizontal  polari-  i 

-V  1 

zation  and  r ,  =  xx  +  yy  or  just  r  evaluated  on  S  .  The  discontinuity  in 

°  1 

D  is  given  by 

J 

1 

AD  =  e^|  'E^(z  =  0‘*')  -  E^(z=  0")| 

1 

] 

(4.7) 

J 

and 

E(z  =  0‘*‘)  =  (1  +Rj^)  (S^) 

(4.8a) 

E(z  «0-)  =  Tj^  l  (S^) 

(4.8b) 

where  Tj^  is  the  Fresnel  (field)  transmission  coefficient  for  horizontal 

polarization.  Combining  (4.8a)  and  (4.8b)  in  (4, 7) yields 

■ 

P 

Ad  =  e  fl  +R.  -  e  T.  1  E,  (S  ) 
oL  h  r  hj  1'  o' 

■*. 

or 

(4.9) 

:l 

Multiplying  (4.7)  by  (4.9)  and  realizing  that  T^  =  1  +  R^  yields 

n  n 

: 

AD-E’  =  -e^,E^^l  +  Rj^)2(G^-l)exp(-j  2k^«rJ  (g.g*) 

1 

(4.10) 

-i 

i 

)  . 

E ' 

Substituting  this  result  in  (4.5)  produces  the  desired  result  for  the 

order  perturbation  field  polarized  in  the  e'  direction 

first 

1 

:i 
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2 

k  ^  E  r  r 


exp(  -  j  2k^Tj^)  dx  dy 


(A. 11) 


where  i'Q'  ■  1, 


cosB-^e  -  sin^  0 

R,  -  - 1_£ - 

"  J - 2~ 

coaO  ■*■  ~  sin  0 

and  so 

„  4(e  -Dcos^e 

(1+R^)^(e^-1)  =  - ^ ^  (4.12) 

l^cos  0  +'^e^  -  sin^O  j 

There  are  several  points  to  note  about  (4.11).  The  derivation  of  (4.11) 
was  considerably  simpler  than  the  Kaylcigh-Rice  approach  [1];  this  is  because 
all  of  the  difficult  work  was  done  in  obtaining  (4.3).  Equation  (4,11)  is  an 
expression  for  the  scattered  first  order  perturbation  field,  a  more  toeaning- 
ful  quantity  than  the  average  scattered  power  when  dealing  with  phase  sensi¬ 
tive  systems.  The  average  of  (4.11)  is  zero  because  <f;  >  =  0j  however,  the 

s 

2 

average  of  <l6*E*e*|  >  or  the  Incoherent  power  is  not  zero*  At  the  beginning 
of  this  section  e'  was  specified  to  be  in  the  same  direction  as  e  ,  thus 
5^E«e*  represents  where  the  double-h  subscript  denotes  horizontal 

polarization  on  transmission  and  reception.  It  is  now  possible  to  examine  the 
consequences  of  cross-polarized  sampling  of  the  scattered  field.  In  this  case 
e*  is  orthogonal  to  e  ;  thus,  e  is  horizontally  polarized  and  e'  is 
vertically  polarized,  l.e. 

S  “  “  sin  4)51  +  cos  4*9 

e*  «  “  -  cos  0  cos  4>x  -  cos  6  sln4)  y  +  sin  0  z 
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Returning  to  (4 >3)  for  this  case.  It  Is  noted  that  AE  Is  still  sero  because 
is  tangential  to  ,  AB  and  AH  are  still  zero  because  there  Is  no 
change  in  across  and  the  conductivity  is  assumed  finite,  and  the  prob¬ 
lem  reduces  to  evaluating  on  .  However,  AD  has  the  direction 

while  E*  is  polarized  in  the  e^  direction,  consequently,  e*e’  *e^*e^*0 
and  there  is  no  depolarization  by  the  surface.  This  is  just  a  confirmation  of 
the  fact  that  first  order  perturbation  theory  does  not  lead  to  a  depolarized 
scattered  field  when  the  rougtmess  is  small  scale. 


4.2.2  Vertical  Polarization 

The  case  of  vertical  polarization  is  a  bit  more  algebraically  involved 
because  AE  is  no  longer  zero  across  ,  Both  magnetic  field  discontinu- 

ities,  AB  and  AH  ,  are  still  zero  for  the  same  reason  as  given  above. 
Thus,  (4.3)  reduces  to 


-  k  exp(  -  j  k 

6‘E*e’  =  — - 


4ii  R  E  e 
o  o 


°L—  J  IAE*D'  -AD»E’]CgdS^ 


(4.13) 


For  vertical  polarization,  it  is  customary  to  use  the  incident  magnetic  field 
Hj  as  the  source.  Thus,  the  incident,  reflected  and  transmitted  magnetic 
fields  on  are  given  by 

Hf  =  ej^exp(  -  j 

H  =  H  R  e,  exp(  -  j  k  ‘r,  ) 
r  ovh 

H.  =  H  T  e.  exp(  -  j  k  ’r,  ) 
t  ovh 


where  R^  and  T^  are  the  Fresnel  (field)  reflection  and  transmission  co¬ 
efficients  for  vertical  polarization.  On  the  corresponding  electric  fields 

are  as  follows; 
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where  q  ■  1,  r  and  and  k^.  specify  the  direction  of  propagation  of  the 

incident,  reflected,  and  transmitted  fields.  Note  that  since  rj.  is  on  , 


which  is  merely  a  restatement  of  the  fact  that  the  angle  of  Incidence  equals 
the  angle  of  reflection  and  Snell’s  law  [9]  is  obeyed.  For  backscatterlng 
and  similar  polarization  sampling  of  the  scattered  field,  the  primed  fields 
are  the  same  as  those  above. 

Although  somewhat  cumbersome  at  this  stage  of  the  development,  it  Is 
desirable  to  split  the  fields  into  components  which  are  tangential  to  and 
normal  to  .  The  reason  for  introducing  this  t  ransformation  is  that  it 
will  be  very  useful  in  the  composite  surface  development  and  it  is  therefore 
benef<  si  to  obtain  some  facility  with  the  technique  on  this  easier  problem. 

The  : ’lal  to  '  is  n  *  z  while  the  tangent  will  be  taken  as  ?  =  n  x  e,  . 

^  n 

Thu.  i>articular  choice  of  T  is  convenient  because  T  •  e  selects  the  com- 

"  V 

ponent  of  e^  that  is  tangent  to  the  surface.  Since  Ae»t  =  0  and  AD»n  =  0 
on  ,  (4.13)  simplifies  to 


\  exp  (  -  j  k^  R) 

47T  ^  .  G 
w  o 


j'  [(AE*n 


(4.14) 


where  e’  has  been  replaced  by  its  equivalent  e^  since  like  polarization 
sampling  of  the  scatte"^  xield  has  been  specified.  The  boundary  conditions 
Ae»t  “  0  and  =  0  should  not  be  discarded  because  they  will  provide  some 

^  A 

useful  relationships.  From  AE»t  =  0  there  results 
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(A.  15) 


[(kjxs^)  +  R^(k_,xa|_)].  X 


(kjXS^)  •? 


while  AS»ft  »  0  yields 


[(kj  xgj_)  +  R^(k^  x?^)]  T„(k^  X  .  a  (4.16) 

where  AE  =  E^+E^-E^  and  AD*Dj+D^-D^  have  been  used  along  with  (4  >15) 

-►  ■> 

and  D  »  \  E  . 

It  is  now  necessary  to  determine  the  field  quantities  inside  the  Integra- 
tion  In  f4.14).  The  quantity  AE*n  can  be  reduced  to  the  following  form 
through  the  use  of  (4.15); 


V: 


Jy  T 

~  H  — ^  (e  -1)  exp  (  -  j  k, ‘r  )  (k  X  e  )»n 
co^r  iXth 


(A. 17) 


while  (4.16)  simplifies  D* ‘n  to 


D-xi  =  -  exp(-j  k^*r^)(fc^xg^).fi  (4.18) 


so 


(AE*n)(D'»n)  =  h  H  T  (c  -  1)  exp  (  -  j  2k  •?  )  ^(k^  *^1  (A. 19) 

oovr  ix'-  cnJ 


Through  similar  manipulations. 


^  (1  -e^)  exp  (  -  j 


(4.20) 


and 


-/r„ 

*  o 


T 


''  exp(- j  k^»rj^)(k^x  ej^)»T  (4.21) 


/r 

r 

42 


(A5-t)(E'-t)  -  (l-c^)exp(-j  [(k^xGj^).^]  (4.22) 


Combining  (4.19)  and  (4.22)  and  completing  the  unit  vector  operations  yields 


(AE‘S)(D'‘n)  -  (AD*t)(E’-t) 


where 


2„2 


V  “V“tr  O] exp (  -  J  2k^Tj_) 


(4.23) 


2c  cos  6 
r 


e  cos  0  +  V  e  -  sin  0 
r  r 


Substituting  this  result  in  (4.14)  and  recognizing  that  E  “  W  ~ 

-  °  f  o  ° 

and  e  "-k.  xe,  ,  the  final  result  is  obtained 
V  i  h  ’ 


k^exp(-J  k  R)  fir~  ,  -1) 

6t.8  .-a °  J^H  — V- 

V  4TrR  ^eov  ^2 

10  e 

r 


[Cr  (Cj.-l)sin^^ j' axp(-j2k^«r^)Cgdxdy 


(4.24) 


Essentially  the  same  remarks  apply  to  the  vertically  polarized  scattered 
field  as  for  the  horizontal  case.  In  addition,  it  should  be  noted  that  if 
(4.11)  and  (4.24)  are  converted  to  o“  or  the  scattering  cross  section  per 


unit  area  according  to 


o“(0,(|>)  “  11m  lira 

R-mjo  a-»« 


J 


where  A  is  the  illuminated  area,  the  result  is  identical  to  the  result 
obtained  by  Peake  using  the  Rayleigh-Rice  approach  [1]. 


■r.v.iik 


4.3  Bistatlc  Scattering  From  A  Dielectric  Surface  With  Small  Roughnase 


For  bistatlc  scattering  the  unit  vectors  specifying  the  directions  of 
Incidence  of  the  unprlmed  and  primed  fields  are  given  by  (see  Figure  4-2) 

-  sin  0^  cos  x  -  sin  6^^  sin  4*^^  9  “  co*  0^  * 

(4.25) 

k!  *  -  sin  6  cos  (b  x  -  sin  0  sin  di  y  -  cos  0  z 
i  s  s  s  ^s  8 

”►  A.  I 

and  k.  *  k  k.  *  k.  *  k  k.  •  The  unit  vectors  specifying  the  directions  of 
ioi  iol  r  j  o 

horizontal  and  vertical  polarizations  for  the  primed  and  unprimed  fields 
are  as  follows; 

e^  =  -  sin(j)^  X  +  coscb^  9 

e,*  =  -  sincb  x  +  cosd)  y 

h  s  s 

(4.26) 

e^  =-cos6j^  coscf)^  X  -  cosO^^  sin<l>^  y  +  sinO^  z 

e*  «  -cosO  cos(b  x  -  cos6  sind)  y  +  sin0  z 

V  s  s  s  s  s 

A  A. 

The  normal  to  the  unperturbed  surface  is  n  =  z  while  the  tangents  to  the 
surface  for  vertical  polarizations  are  given  by 

T  =  n  X  a.  t'  =  n  X  e  ’  (4.27) 

n  n 

Equations  (4.25)-(4.27)  d, scribe  the  basic  quantities  that  will  be  required  in 
this  section. 


4.3.1  Horizontal  Polarization 

Since  AB  and  AH  are  zero,  because  there  is  no  change  in  across 
and  the  conductivity  is  finite,  (4.3)  becomes 


6lE 


f  - 


exp(-Jk  R)  f  ^ 

o  o  J 


•E'  ]  d  S 
s  o 


(4.28) 
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I 

I 

( 

t 

] 

i 

f 


Furthermore,  AE  is  also  zero  across  S  because  E  is  tangential  to  . 

o  o 

Thus ,  (4,28)  becomes 


n 


exp  (  -  J  R) 

AttR  E  e 
o  o 


jld- 


E'  1  C,  iS^ 


(4.29) 


On  the  surface  S  ,  the  incident,  reflected,  and  transmitted  electric  fields 
o  ’  ’ 

(unprimed  and  primed)  are  as  follows; 


^1  = 

E^  exp(- jk^.r^)a,^ 

“ 

Eoexp(-jit;.?^)S^' 

-► 

->  .  J. 

. 

E  * 
r 

F,  * 
r 

e 

ft  .  ft  , 

->■ 

.  V  Vk. 

E  T,  exp(  -  j  k  ‘r  )e. 
oh  ‘  ^  t  X  e 

ft  «  ft  . 

EoT^^'exp(-jkJ-r^. 

where  also  on  the  surface  S  , 

o 

(4.30) 

-►  -> 

k .T  . 

=  k  ‘r  =  k  ‘r 

, 

•r, 

■=  k  •  r  ■  k  '  •  r 

(4.31) 

IX  r  _L  t  X 

i 

JL 

r  JL  t  X 

and 

the  same 

:  notation  as  introduced  earlier 

has 

been  continued.  The  fields 

-y 

AU 

and  E* 

on  S  are  as  follows; 
o  ’ 

Ad  =  e  E  (1  +  R  -  e  T,  )  exp  (  -  j  k .  •  r  )  S. 
oo  h  rh  *^  ''iX  h 

E’  =  E^(l  +R^‘  )  exp  (  -  j  k’^4^) 


Using  the  fact  that  1  +  =  Tj^  ,  the  product  Ad»E'  becomes 

AD•■^^  =  E^^  (1  +  Rj^)  (1  )  (1  -  c^)  )  exp  [  -  j  (k^  +1^^'  )  j 

and  substituting  this  result  in  (4.29)  yields 


(.S^E*g,' 


exp(-jk 


AttR 


- —  J J'exp[-j(k^+kp.rJCg  dxdy 


(4.32) 
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This  result  csn  be  easily  translated  into  the  angles  (0  ,(j>  )  and  (0  ) 

IX  S  H 

by  the  use  of  (4 >25)  and 


®h  ®h 


cos (A.  -  A  ) 
s 


along  with 


2  cos  6 , 


2cos  6 


cos  0.+ -  sin^  6. 
i  r  J 


1  +  R,* 

h 


0  +  -  sln^  0 


cos 


For  ease  of  comparison,  It  should  be  noted  that  most  results  similar  to  (4.32) 

A 

express  the  direction  of  scattering  as  k  which,  in  the  above  notation,  is 

s 

-  k^’  .  Equation  (4,32)  yields  .  For  cross  polarized  sampling  of  the 

scattered  field,  the  problem  becomes  somewhat  more  Involved  and  it  will  be 
discussed  in  Section  4.3.3.  It  ihould  be  noted  that  when  <1)^  =  (j)^  (backscat- 
tering) ,  (4.32)  reduces  to  the  result  obtained  in  Section  4.2,1. 


4.3.2  Vertical  Polarization 

Since  n«E  and  t*D  are  discontinuous  across  the  unperturbed  surface, 
(4.3)  reduces  to 


k^  exp(-jk^R) 

4'itRE  e 
o  o 


J*  |^(AE»n)  (o’*  a)  - 


(AD»t)  (E’  •!'  )  (l  •l' 


i;  d  s 

s  o 


(4.33) 


A  A  -  A 

where  the  scalar  product  must  be  included  because  the  unit  vectors  x 

and  T*  are  not  necessarily  parallel,  e.g,  see  (4.27).  The  Incident,  reflec¬ 
ted,  and  transmitted  unprimed  and  primed  magnetic  field  quantities  are  given  by 


«i 


H^exp(-jkj-r_^)ej; 


H  ’ 

r 


St 


H  R  '  exp( -J '•r  >e  * 
ov  *^  ''r  ih 


=  H^T;exp(-jit^‘.r^)e^' 


(4.34) 
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while  the  corresponding  electric  fields  are 


V^^o  - 

—  k  X  H 
e  i  i 


-  k^’  X  H^' 


E  =  -  k  X  H 

r  1C  r  r 
o 


->  J  ‘'o  -  -> 

E  =  -  W— ^  k  X  H 

t  1  c  c  t  t 


"t’  = 


—  k  '  X  H  ' 
e  r 

o 


V  "  "t 


(4.35) 


e  e 
o  r 


and  D  =  e  E  ,  D'  =e  E*  .  The  same  notation  introduced  in  Section  4.3.2  is 
continued  here.  Frv^m  the  boundary  condition  that  the  tangential  component  of 
the  electric  field  be  continuous  across  ,  i.e.  AE*t  =  0  and  6E**T*  =  0  , 


the  following  relationships  result; 


T 

^  V 


(4.36) 


^r 


'S'  xg/).9’ 


Similarly  from  the  continuity  of  the  normal  component  of  the  D-field  across 

S  ,  i.e.  AD*n  =  0  and  AD'»n  =  0  ,  there  results 
o 

(k.  xS  ).fi  +  R  (k  xg  ).fi  =  y'T'  T  (k  xg  )»n 
'i  li  v  r  u  rvt  h 


(4.37) 


(k;  xg  •).g  +  R  •  (k  •  xg  *).fi  =  »^  T  ’  (C  '  xg  •).g 
i  h  v  r  h  rvt  h 


Using  (4.37)  to  simplify  the  expressions  for  AE»fi  and  D' ‘n  yields 


(4.38) 


Al*fi 


-VF 


e”  Ho  ^  -  l)exp(  -  .1  -T^)  (kj.  X  .Cl 

'^r 


and 


D'.fi  -  -  ^  e  H  v-e  T  •  exp(  -  j  kj'r,  )(k  '  xg  *).fi  (4.39) 

ooo  rv  h 

In  a  similar  fashion,  (4.36)  Is  used  to  simplify  the  expressions  for  AD*t 


^  t 

and  E  ‘T  with  the  following  result: 


AD*T 


Ho  — exp(-ikj‘rj^)(kj^x|^).T 

r 


(4.40) 


and 


E’«t'  = 


-VF 


T  *  ^ 

Ho  —  exp  (  -  j  •T)  (k  J  X  )  .^ • 

}/c~ 

r 


(4.41) 


Substituting  (4.38)  -  (4.41)  in  (4.33),  noting  that  t»t'  =  cos(<l).  -(I)  )  , 

X  s 

and  simplifying  the  unit  vector  operations  yields  the  following  result  for 


•  ; 
V 


^  H 

▼  e  0 


exp (-jkoR) 


6*E*e  '  -  -  ^  .. 

V  *  e  o  4ttr 

o 


’  T  '  — —  /e  slnO,  sin0  +  /(c  -sin^O . )  (e  -sin^G  ) 
vv^Z^^ri  s  "r  ir  s 


cos (4)^  -  (j)^)^  J* J' exp  [-  j  (k^  +k^’  )  dxdy 


(4.42) 


where  E  “  H  /p  /e  and 
0  ooo 


T  -  1  +  R  = 

V  V 


2  e  cos6 . 
r  i 


e  COS0,  +  -  sln^0. 

r  I  r  1 
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-A:.  Li,.  .  ;  .,1 


.awwanWWIKAWi 


and  it  is  easily  shown 


and  T  '  =  T  ( 0  ->-0  )  .  liquation  (4.42)  is  fi'li  , 

V  v'  1  s  ^  w* 

to  reduce  to  the  result  for  backscattcring  in  Section  4.2.2.  when  '!> .  *  <t* 

X  8 

and  0=6  .  Once  again  it  should  be  emphasized  that  the  most  difficult 

1  s 

part  of  obtaining  (4.42)  is  evaluating  the  terms  (k^xe^^)*!  and  (kJxg^*)»T' 
in  (4.40)  and  (4.41).  As  noted  previously  this  is  a  consequence  of  the  fact 
that  the  difficult  analysis  was  finished  once  (4.3)  was  derived  and  the  actual 
evaluation  of  (4.3)  is  very  straightforward.  Finally,  comparing  the  a® 
values  resulting  from  (4.32)  and  (4.42)  with  the  corresponding  results  obtained 
from  the  Rayleigh-Rice  theory  [10]  shows  complete  agreement. 

4.3.3  Cross  Polarization 

As  shown  in  Sections  4.2.1  and  4.2.2,  depolarization  for  scattering  by 
small  scale  roughness  is  a  second  order  effect  ^  the  plane  of  incidence. 

If,  however,  the  scattered  field  outside  of  the  plane  of  incidence  is  computed, 
it  will  be  found  to  have  a  nonzero  cross  polarized  component.  This  result  is 
simply  a  consequence  of  the  fact  that  the  unit  vectors  and  e^'  are 

not  fixed  with  respect  to  the  surface-centered  coordinate  system  and  they 
change  their  directions  as  the  observation  point  moves  out  of  the  plane  of 
Incidence.  This,  of  course,  is  a  purely  geometrical  effect  and  it  has  nothing 
to  do  with  any  change  in  tlie  basic  scattering  mechanism. 

The  derivation  of  the  results  follows  essentially  the  same  pattern  as 
set  forth  in  the  previous  sections.  There  is  one  point  that  should  be  noted 
because  it  simplifies  tlie  algebra  somewhat.  For  the  case  of  the  incident 
field  liorizon tally  polarized  (Cj^)  and  the  scattered  field  vertically  polarized 
(@^' )  ,  the  unprimed  field  quantities  should  be  obtained  from 

exp(-j  k^»r^)e|^  while  the  primed  fields  should  be  obtained  from 
=  H^exp(-j  the  incident  field  vertically  polarized  (e^) 
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and  the  scattered  field  horizontally  polarized  (^j^' )  •  the  unpriraed  fields  are 
obtained  from  exp(  - j k^*r^)ej^  while  the  primed  fields  are  to  be  de¬ 
rived  from  ®  E^exp(-j  •  This  approach  is  consistent  with  the 

technique  of  obtaining  all  field  quantities  from  the  horizontal  (e,  or  e  ' ) 

n  h 

field  for  planar  surface  reflection. 

With  E  ®  H  /y  /£  ,  the  following  expressions  for  6*E,  ,  and  6* E  ,  , 

O  O  O  O  or 

result ; 


.  k  exp(-jk  R) 

,  -  s'i-r  -  -g-  ■-  ■?-  E 

nv  V  ttR  0 


cos9.  cos6  sln(d>  -<!).)(£  -1) 
i  s  '^s  i  r  ' 


-dvT^ 


sin^0 
r  s 


(cos0^  +k4 


'  +sin^0 .  )  ( e  cos  0  +  v4  +sin^0  ) 
r  1  r  s  r  s 


(4.43) 


k^exp (-jk^R) 

Tr 


^  - > 

cos0,cos0  sin((i)  -(J>, )  (e  -1) 

X  s  ^s  i  r 

<  - ; —  - -  - 

(cos0^+V^^-sin^O^)  (e ^cos6 £f  ^-s in^ 0^ )J 


^^e^-sin  0. 


where,  in  summary,  the  angles  are  defined  in  Figure  4.2,  r^^  *  xx  +  yy  ,  and 

and  are  defined  as  k^k^  and  *  respectively,  where  k^  and 

A 

kj*  are  given  in  (4.25). 

When  comparing  (4.43)  and  (4.44)  with  the  cross  polarized  scattered  fields 
resulting  from  the  Raylelgh-Rlce  approach  [10],  (t>^  should  be  set  equal  to  7T  . 
The  expression  for  6*E^j^,  agrees  with  the  results  in  [10,  pg.  706].  The 
expression  for  6*E^^,  is,  however,  the  negative  of  the  coefficient  in 

(10,  pg.  706,  eqn.  9.1-69],  Normally  this  difference  is  not  important  because 
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is  squared  and  then  averaged  to  find  the  incoherent  power.  However, 

if  one  is  dealing  with  circular  polarization  the  sign  does  become  critical. 

To  resolve  this  issue,  a  special  case  can  be  constructed  whereby  should 

agree  with  a  or  6* E,  ,  =  6* E  .  This  special  case  involves  taking 
vv  hv  vv 

<b  -  Tf  and  d)  =  0  and  0  =  0  in  tho  expression  for  a  or  6^E  ,  and 

1  s  S  w  w 

comparing  this  result  with  and  for  41^  =  it  ,  <|>g  =  Tr/2  ,  and 

0=0.  In  this  special  case  both  6^ E  .  and  E.  1  should  be  polarized 

in  the  -  x-dlrection  .  Comparing  (4.42)  and  (4.43)  for  this  special  case  shows 

that  Indeed  6'e  i  =  6^ E.  ,  .  However,  evaluating  a  and  a,  from  [10] 

vv'  hv'  ’  ®  w  hv  ^ 

results  in  ;  consequently,  there  does  appear  to  be  a  sign  error 

in  the  expression  for  and  (4.43)  is  correct. 

This  section  completes  the  development  for  scattering  from  a  dielectric 
surface  having  only  a  small  scale  roughness.  Once  again  it  should  be  emphasized 
that  the  purposes  of  Sections  4.2  and  4.3  are  (1)  to  check  the  Burrows  pertur¬ 
bation  approach  against  the  conventional  Rayleigh-Rice  results  and  (2)  to  il¬ 
lustrate  the  actual  mechanics  of  evaluating  the  Burrows  expression  for  the 
first  order  perturbation  field.  Hopefully,  this  latter  purpose,  if  achieved, 
should  considerably  simplify  the  transition  to  the  composite  surface  case. 

4.4  Bistatic  Scattering  From  A  Dielectric  Surface  With  Composite  Roughness 
For  small  scale  roughness  superposed  on  a  planar  surface,  the  Burrows 
perturbation  formula  (4.3)  is  particularly  easy  to  evaluate.  This  results  from 
the  fact  that  one  deals  with  the  fields  on  an  infinite  planar  surface  and,  for 
such  a  surface,  the  fields  are  easily  described  and  related  through  the  Fresnel 
coefficients,  Snell's  law,  and  the  equality  of  the  angles  of  incidence  and 
reflection.  For  a  composite  surface,  the  unperturbed  surface  is  not  planar 
but  it  is  assumed  to  be  very  gently  undulating.  More  specifically,  the  unper¬ 
turbed  surface  is  actually  defined  such  that  it  contains  no  spatial  frequency 
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components  vhlch  are  smaller  than  Yk^  ,  where  y  is  a  constant  wlilch  Is 
greater  than  unity.  Of  course,  it  is  desirable  to  have  y  as  large  as  pos¬ 
sible  but  this  is  not  always  practical  since  the  small  scale  height  must  satisfy 
4k^  <<  1  [6]  .  However,  if  y  can  be  made  sufficiently  large  then  che 

scattering  from  the  unperturbed  surface  can  be  treated  using  physical  optics. 
Physic?  nf'lcs  assumes  that  the  surface  may  be  considered  to  be  locally  planar 
and  ti;t  fields  on  the  surface  can  be  accurately  approximated  using  Fresnel 
th>..'ry.  This  approach  is  recognized  to  be  essentially  the  same  as  the  small 
scale  roughness  on  a  planar  surface  problem.  The  one  Important  difference  is 
that  for  the  gently  undulating  unperturbed  surface,  the  local  normal  is  no 
longer  entirely  z-directed  and,  in  fact,  depends  upon  the  slopes  of  the  large 
scale  surface.  This  means  that  one  must  construct  a  local  coordinate  system 
on  the  undulating  unperturbed  surface  and  compute  the  surface  fields  required 
in  (4.3)  in  terms  of  this  system.  This  must  be  done  for  both  the  unprimed  and 
primed  fields  because  they  have  different  angles  and  directions  of  incidence 
for  the  general  bistatic  case. 

For  the  unprimed  fields,  the  important  unit  vectors  are  and  n^^ 

which  is  the  normal  to  the  large  scale  or  unperturbed  surface.  These  two  quan¬ 
tities  are  important  because  they  form  the  local  plane  of  incidence.  One  next 
constructs  unit  vectors  and  e^^  which  are  orthogonal  and  parallel, 

respectively,  to  k^^  and  n^^  .  These  unit  vectors  are  also  horizontally  and 
vertically  polarized,  respectively,  with  respect  to  the  local  plane  of  inci¬ 
dence.  Any  arbitrarily  polarized  unprimed  incident  field  can  now  be  decomposed 
into  components  parallel  to  e^^  and  since  the  incident  field  must  be 

transverse  to  k^  .  The  unprimed  field  quantities  required  in  (4.3)  can  then 
be  computed  as  in  the  previous  sections.  The  exact  same  construction  of  e^^j^ 
and  and  the  decomposition  of  the  primed  incident  field  must  be  performed 
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in  order  to  compute  the  primed  fields  required  in  (4.3).  Fortunately,  this 
is  easily  accomplished  by  simply  changing  6.  to  6  and  <{).  to  (j)  in 
the  unprlmed  quantities. 

Before  getting  into  the  actual  details,  there  are  a  few  other  points  that 
should  be  noted.  All  of  the  above  noted  manipulations  are  going  to  lead  to 
the  following  changes  in  6^E  obtained  in  equation  (16)  of  [6].  First,  the 
factor  is  going  to  depend  on  ,  the  angles  and  directions  of  inci¬ 
dence  and  scattering  (0^  ’  *^i  ’  *  ^s^  *  slop®®  large  scale  or 

unperturbed  surface  only  other  change  is  that  the  exponen- 

tial  inside  the  surface  integral  will  become  expt  -  j(k^ +k^’ ) 'rj^l  where 
=  XX  +  yy  +  z  because  of  the  generalization  to  blstatic  scattering. 
Except  for  correcting  [6]  to  properly  Include  shadowing,  as  detailed  in  Section 
2,  all  other  aspects  of  the  solution  presented  in  [6]  remain  the  same.  Com- 
bining  6*E  from  this  analysis  with  Sancer's  result  [11]  for  essentially  6  E 
yields  the  total  scattered  field.  Furthermore,  it  should  be  expected  that  for 
backscatterlng  the  dielectric  nature  of  the  surface  should  have  an  almost 
negligible  effect  upon  the  wavenumber  at  which  the  surface  height  spectrum 
is  partitioned  into  large  and  small  scale  sub-spectra.  Finally,  because  there 
is  no  discontinuity  in  magnetic  properties  across  the  unperturbed  surface  and 
the  conductivity  is  assumed  to  be  finite,  Ab  and  Ah  in  (4.3)  will  be  zero. 


This  fact  holds  true  regardless  of  any  tilting  of  the  locally  planar  surface. 
Tlie  first  task  at  hand  is  to  constr  ict  e^^^  and  .  Since  e^^^^  is 

A  ^ 

orthogonal  to  both  k^  and  ,  it  is  given  by 


(4.45) 


where  k.  =  k,  x  +  k.  y  +  k^  z  and  the  k,  ,  q  =  x,y  and  z  ,  are  obvious  from 
j.  lx  ly  iz  iq  ^  ’ 

(4.25).  The  normal  to  the  large  scale  surface  is  given  by 
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(4.46) 


>/l  +  4  ^  ^ 

Jlx  8,y 


n.  X  +  n.  9  +  n-  z 
!.x  iy  •'  Iz 


Expanding  the  cross  product  In  (4.45)  yields 


=  ^x  "  ^  ‘'£y  y  ^  '‘£x 


(4.47) 


where 


(k.  n-  -  k,  n.  ) 
iy  X.Z  iz  £y 


l^l  X  SJ 


(k.  n.  -  k.  n.  ) 
Iz  £x  ix  £z 


1^1  X  (ijl 


(k.  no  -  k,  n.  ) 
ix  iy  iy  ix 


l^i  X  fiji 


and  e.v  is  completely  determined.  For  the  unit  vector  e-  ,  tl\e  expres- 


5  *1 

^  4 


slons  are  more  involved  because  must  be  in  the  plane  form  by  k^  and  n^^ 
and  also  orthogonal  to  e^^j^  .  For  the  retlected  and  transmitted  fields  k^ 

A  A  ^ 

goes  to  k  and  k  «  This  will  not  change  the  direction  of  e.,  because 
r  t  x»n 

k^  ,  and  are  all  coplanar.  Tills  will,  however,  alter  the  direction  of 
e^^  ;  this  is  easily  understood  by  noting  that  k^  ,  e^j^  ,  and  form  a 

mutually  orthogonal  triad  of  unit  vectors.  Thus,  if  does  not  change 

A 

direction  but  k^  does,  then  e^^  must  necessarily  change  direction.  What 

this  means  is  that  we  must  find  a  new  e^  for  each  value  of  k  .  This  is 

iv  q 

easily  done  by  the  following  equality; 


iv  ih  q 


(4.48) 


where  q  -  i,  r,  and  t  .  Note  that  it  is  not  necessary  to  divide  the  rhs  of 
(4.48)  by  the  magnitude  because  it  is  unity,  i.e.  and  Ic^  are  mutually 

orthogonal  by  (4.45).  The  unit  vector  e-'*  may  also  be  written  as 


(4.49) 
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where  6^^.  -  /E^  for  q  -  i  and  r,  and  -  1  .  Expanding  (4.52)  yields 


*1  ■  (-  <  'tv  i  'ti, ) 

St  •  V¥^  (-  <  'ti '  'th) 


From  Fresnel  theory,  the  Cj^j^-component  of  is  equal  to  R  times  the 

component  of  »  so 


E*^  =  R  E* 
V  V 


(4.53) 


Similarly,  the  component  of  is  equal  to  times  the  e^^j^-component 


of  ,  so 


E^ 


(4.54) 


Equations  (4«53)  and  (4.54)  can  now  be  used  in  (4.50)  to  express  all  of  the 
fields  in  terms  of  E^  and  ,  i.e. 

E  =  E?;  e  +  E^  e/ 

1  h  ih  V  xv 


E  =  R.  E^  e,.  +  R  E  e„ 
r  hjj^  h  ih 


(4.55) 


E..  =  T.  E,  e-,  +  „1  t 

t  h-  h  ih - E  On 

%  V  Jlv 

r 

«The  subscripts  on  R^,  Rj^,  T^,  Tj,  means  that  the  angles  in  the  appropriate 
Fresnel  formulas  must  be  defined  with  respect  to  the  normal  to  the  large  scale 
surface,  h 
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i  i 

where  Ej^  and  are  given  by  (4*52)  with  q«i  .  The  corresponding  D-fields 

are;  D.  =  c  K,  ,  D  •=■  c  E  ,  «  e  c  .  The  printed  fields  may  be  obtained 

from  (A  55)  and  (4.51)  by  replacing  by  k^'  ,  by  e^j^.  by  and 

by  changing  e  to  whatever  scattered  field  polarization  is  to  be  sampled, 
a 

r^ay,  . 

The  appropriate  form  of  (4.3)  la 


b 


k^2exp(  -  ,1  k^R) 

AttR  E  e 
o  o 


S 

o 


^  d  S 
8  o 


(4.56) 


where  the.  shadowing  factor  has  been  temporarily  omitted  from  the  integrand 

Kt.r- 


since  it  can  be  added  at  the  end  of  the  development.  The  (AE'fij^)  is  equal 
to  (E^ +E^ -E^)»nj^  or  using  (4.55)  and  noting  that  =  0  , 


Equation  (4.57)  can  be  simplified  somewhat  by  using  the  relationship  that 

•>  A 

results  from  AU*n  =  0  .  The  final  result  is 


1 

AE*fi.  -  E  — - -  T 


For  0’»n 


a  ’ 


i^*h'  =  E^'  e  i/r~  T'  (o.y*n.) 

a  V  or  v^  Kv  * 


(4.58) 


(4.59) 


so  the  product  of  (4.58)  and  (4.59)  can  be  written  as  follovjs 

(S'  -Bj)  -  [.J  ]  Cp (Ep-1)  T'  T  (Sj^'  -Sj)  (4.60) 
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It  can  be  shown  that 


e^-a 
Iv  A 


^  t'  .A 

®J!,v  ”)l 


where 


so  (4.60)  becomes 


sin  9, 


'l-(-k^«nj^)' 


5in0,'i=  Vl-(-k^*nj^) 


.xp[. 

(4.61) 

For  the  temalnlng  tern  in  (4.56),  the  Important  parta  are  (AD)^  and 
(!■)  where  the  p-aubacript  denotea  tangential  to  the  large  scale  aurface, 

l.e.^the  normal  component  o£  D  la  contlnuoua  acrosa  the  boundary  so  AD-Oj-O 

The  tangential  components  o£  the  A5-Eleld  can  be  Eound  by  decomposing  A?  Into 
components  directed  along  6^^  and  1,^  .  where  »  Sg  ■<  .  l.e. 


ilsing 


(AD)p  =  +  (AD-Tj^)Tj^ 


(4.62) 
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and  simplifying  this  expression  with  the  aid  of 


it  t  A 


which  results  from  (AE)^  =  0  yields 


(4.63) 


The  appropriate  expression  for  (^*)p 


p  "h  h^  Ah 


''  ^ 


T,;  +  E^'  <^V 


(4.64) 


since  on  the  unperturbed  surface  E*  =  E^+E^  *  E^  .  The  relationship  for 

(^)  (E')  is  obtained  by  taking  the  dot  product  of  (4.63)  with  (4.64).>  Com* 
P  P 

bluing  this  result  with  (4.61)  and  substituting  into  (4.56)  yields  the  follow¬ 
ing  result; 


(5‘e 


E  k 
o  o 


exp(-jk^R) 


ab 


TTR 


// 


^ab^^Ax’^Ay 


)I(x,y)exp  1^- J  (kj.+kj^' )  ‘r^  ] 


(4.65) 


where  6'e  .  is  the  scattered  first  order  perturbation  field  for  an  incident 
ab 

polarization  e  and  scattered  polarization  e.  ,  I(x,y)  is  unity  on  the 
illuminated  parts  of  the  surface  and  zero  for  the  shadowed  parts,  and 


(e^-1) 


X  y 


rfA  X  V  >A  ^  1.  V 

(e  *6.  )(e, ‘e.  ) 

a  Jlv' '  b  8.V 


T  T  ' 
'*1 

2 


V  e. 


sin  0^^  sin  0^^ 


(AAt  A  A 

sin  sin 


k  •k/  -  (np*k.)(no*k/)  ^  ^ 

sinGj^^  ainej^^ 


/br-sin 


-A  ^  i*  V  /A  ^  \  m  (  ^ 


sin  9j^j  sin  O^'j 


(4.66) 


For  convenience,  the  above  terms  are  sutianarized  below 


^  i  A 

S  -  Polarization  of  the  incident  electric  field  (E,  = 
a  X  o  a 


e^  ■  Polarization  of  the  scattered  electric  field  (E^  “ 


A  A  -  A 

liixjLifiii:!:! 
^  hi  *  hy 


‘‘i  ■  ‘^oh 


-  k  k.' 
i  o  i 


f'  ^  n 

k^  »  -  sin  cos  ({>j^  x  -  sin  G^  sin  <}>j^  y  ”  cos  0^  z 

k,'  «  -  sin  G  cos  <J)  x  -  sin  0  sin  (|>  y  ”  cos  0  z 
i  s  s  s  s '  s 


A 

A  A 

i  *  ^  t 

k’  xn^ 

®ilh  " 

iT  1 

IH  "i' 

sin  G 

1  A  A  • 

U  “ 

aln  0^^  - 

-> 

^S,  " 

X  X  +  y  y  +  r,^  Z 

nos  Oy  . 

'1,  ‘  ®ih  “  ''i 


I-Cjxnj 


(-kj)-Su 
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5.0  A  USEFUL  RELATIONSHIP  FOR  THE  JOINT  SLOPE  PROBABILITY  DENSITY  FUNCTION 
5 . 1  Background 

The  Incoherent  power  scattered  in  and  about  the  specular  direction  de¬ 
pends  upon  the  Joint  probability  density  function  (jpdf)  of  the  large  scale 
slopes.  The  jpdf  for  the  large  scale  slopes  is  also  important  in  determin¬ 
ing  the  degree  of  "tilt"  or  k-space  broadening  imparted  to  the  small  scale 
Bragg  scatterers.  Ideally,  one  would  like  to  measure  the  jpdf  for  the  large 
scale  slopes,  the  roughness  spectruia  of  the  small  scale  heights,  and  the  com¬ 
plex  dielectric  constant  of  the  surface  in  order  to  predict  the  average  scat¬ 
tering  properties  of  a  specified  section  of  terrain.  That  is,  these  surface 
measurements  would  be  substituted  in  the  rough  surface  scattering  model  which, 
in  turn,  would  provide  an  estimate  of  the  average  coherent  and  incoherent 
scattered  power.  From  a  practical  point  of  view,  measurements  of  the  jpdf 
for  the  slopes  and  the  small  scale  roughness  spectrum  are  very  difficult  to 
obtain  and  the  difficulty  increases  as  the  radar  or  electromagnetic  wave¬ 
length  decreases.  For  example,  in  the  case  of  an  L-band  system  with 

=  30  cm,  the  small  scale  part  of  the  scattering  model  will  require  sur¬ 
face  height  spectral  measurements  of  surface  undulations  having  wavelengths 
of  less  than  about  90  cm  because  ~  /2)csc  0.  for  backscatter. 

For  the  large  scale  features  of  the  surface,  the  jpdf  for  the  slopes  repre¬ 
senting  surface  features  having  spatial  wavelengths  greater  than  about  90  cm 
is  required.  Obviously,  spectral  information  on  the  small  scale  features  is 
going  to  be  the  most  difficult  to  obtain.  However,  even  the  jpdf  for  the 
large  scale  slopes  is  going  to  be  difficult  to  estimate.  It  is  not  unrea¬ 
sonable  to  expect  that  we  can  obtain  measurements  of  the  jpdf  for  the  large 
scale  heights  and  even  the  correlation  function  for  the  heights  at  least  some¬ 
what  close  to  the  90  cm  spatial  resolution.  However,  this  information  must 
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somehow  be  translated  Into  the  jpdf  for  the  surface  slopes  and  this  Is  where 


the  difficulty  comes  in.  We  Ignore  the  possibility  of  a  direct  measurement 
of  the  jpdf  of  the  large  scale  slopes  for  arbitrary  terrain  because  such  a 
task  appears  to  be  too  difficult  to  even  contemplate. 

The  question  basically  bolls  down  to  the  feasibility  of  translating  or 
converting  measurements  of  the  Jpdf  for  the  surface  heights  into  the  jpdf  for 
the  surface  slopes.  The  purpose  of  this  section  is  to  point  out  an  analyti¬ 
cal  means  for  accomplishing  this  transformation  and  suggest  that  the  scheme 
be  attempted  on  an  experimental  basis.  The  relationship  is  not  new  and,  in 
fact,  results  from  some  earlier  rough  surface  scattering  analysis.  However, 
it  has  apparently  gone  unnoticed  at  least  insofar  as  it  applies  to  this  very 
real  world  problem  of  translating  the  height  jpdf  into  the  slope  jpdf. 

5.2  The  Transformation 

Perhaps  the  oldest  approach  to  estimating  the  quasl-specular  incoherent 
power  scattered  by  a  rough  surface  is  now  called  the  autocorrelation  approach 
Basically,  one  assumes  the  validity  of  physical  optics,  interchanges  the 
order  of  spatial  Integration  and  ensemble  averaging  in  the  expression  for  the 
scattered  power,  and  assumes  Gaussian  surface  statistics  with  the  final  re¬ 
sult  that  the  average  scattered  power  is  dependent  upon  the  behavior  of  the 
surface  height  correlation  function  near  |Ar|  =  0[1]  .  In  the  mid-60's, 
Kodls  [2]  showed  that  the  average  scattered  power  could  alternatively  be 
interpreted  in  terms  of  the  number  of  specular  points  on  the  surface  and  the 
absolute  radii  of  curvature  at  the  specular  points.  Barrick  [3]  subsequently 
linked  these  two  approaches  in  the  high  frequency  limit  where  both  are  valid. 

In  the  process  of  establishing  the  similarity  between  the  autocorrela¬ 
tion  and  specular  point  approaches,  Barrick  obtained  a  relationship  between 
the  jpdf’s  of  the  surface  heights  and  the  st5rface  slopes.  In  particular, 
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C  the  x  and  y  surface  slope  components  and 

X  y 

the  joint  characteristic  function  for  the  surface  heights 


Cl(Xl.yi)  and  ^2^*2 ’^2^ 


^x^y'  ^  ^z  '  4vi^  k  ^{^JJ  ^1 

Q  ^  —00 


y.  (k  =k  q  ,k  «*-k  q  ;A  ,A  ) 
X  o  z  y  o^z  X  y 


•  exp(j  ko^jj  Ax  +  j  k^q^  Ay)dAx  dAy 


(5.1) 


where  Ax  =  x^  -  X2  and  Ay  -  y^^  -  y^  .  In  (5.1)  the  quantities  aad 

q^  are  limited  as  follows;  |q^|  <1  ,  i  =  x,y,z  .  It  should  be  noted  that 
$_  (•)  is  an  implicit  function  of  the  surface  height  correlation  function; 

this  is  how  the  <Ax,Ay)  variation  comes  about  in  (5.1). 

If  an  analytical  form  for  _  is  available  then  (5.1)  can  be  used 

C1C2 

directly  to  obtain  the  Jpdf  for  the  slopes.  In  cases  where  $_  „  (•)  is 

^1^2 

obtained  from  measured  data,  It  is  not  Immediately  obvious  that  (5.1)  is  of 
any  practical  use  since  the  behavior  of  the  joint  height  characteristic  func¬ 
tion  will  not  be  known  in  the  limit  of  k  and  k  .  However,  consider 

X  y 

the  following  reasoning  as  a  means  for  obtaining  estimates  of  P  _(*,*). 

^x> 

Since  $  is  the  two-dimensional  Fourier  transform  of  the  jpdf  for 

the  height,  it  can  be  obtained  numerically  by  using  a  Fast  Fourier  Transform 
(FFT)  on  the  measured  jpdf  height  data.  The  result  of  this  operation  will  be 
denoted  by  ^^^^^2  *  measurement  noise  and  particularly  quantlza- 

A 

tion  noise  in  the  measured  height  jpdf  data,  $  will  be  limited  to  values 

44 

less  than,  say,  k  <  K  and  k  <  K  .  The  maximum  value  of  k  that  can 
X  —  X  y  ~  y  o 

be  achieved  in  (5.1)  is  therefore  max(K  /q  ,K  /q  )  .  If  q  is  small  then 

X  z  y  2  z 

the  resulting  maximum  value  of  k  can  be  very  large.  The  transform  variables 


in  (5.1)  will  be  given  by  K  q  /q  and  K  q  /q  which  may  also  be  large. 

XX  y  y  ^z 


depending  upon  q  and  q 
X  y 

puted  using  the  following; 


Thus,  if  q  is  near  zero  and  P_  ..  is  com- 
z  C  C 

X  y 


?  5 
X  y 


»K  ,  k  -K  ;  Ax,, Ay) 


-L 
K  q 


/  K  q  K  q  \ 

Ax  +  j  ^  Ay  j  dAx 


dAy 


(5.2) 


^  is  a  sufficiently  good  estimate  of  P  as  to 
^x^’y  ^x^y 


it  may  turn  out  that  P^ 
be  useful  in  the  scattering  model.  Unfortunately,  this  approach  breaks  down 


when  q  or  q  =  0  ;  however,  it  may  be  possible  to  get  close  enough  to 
X  y 

q  =0  or  q  =0  to  infer  the  behavior  of  P  along  these  lines  in  the 

^  ^  ’xS 

q,q  “plane.  The  limits  on  the  integrations  in  (5.2)  symbolically  denoted 
X  y 

as  ±  L  ,  will  be  determined  by  the  correlation  length  of  the  surface,  i.e, 
the  separation  distance  for  which  the  surface  height  correlation  function 
is  essentially  zero. 

An  alternate  approach  to  estimating  P  is  to  examine  the  asymptotic 


behavior  of  ‘t> 

possible  to  generate  an  asymptotic  functional  dependence  of  $ 
and  k 


„  „  as  k  -^K  and  k  ->-K  .  From  this  behavior,  it  may  be 
^1^2  x  X  y  y  *  ^ 

r  r  O" 

By  repeating  this  procedure  for  different  values  of  Ax  and  Ay  , 


it  might  be  possible  to  also  generate  or  build-in  the  functional  dependence 

^  A 

of  <I)  on  Ax  and  Ay  .  In  this  manner,  the  dependence  of  0  upon 

k  ,k  ,A  ,  and  A  is  obtained  at  least  in  the  limit  of  moderately  large  k 
X  y  X  y  X 

and  ky  .  This  functional  form  cotild  then  be  transformed  according  to  (5.2). 
The  major  problem  her*  is  that  the  accuracy  of  the  result  will  depend  directly 
upon  how  precisely  the  surface  height  correlation  function  is  known  near 
Ax = 0  and  Ay = 0  .  This  statement  results  from  the  fact  that  the  behavior 
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of  the  v.ransfonn  of  a  function  as  k-*-'®  Is  directly  determined  by  the  behavior 
of  the  function  as  Ar-^-O  [4]. 


The  problem  of  converting  height  jpdf  data  into  slope  jpdf  results  is 
definitely  not  easy.  Even  with  the  use  of  (5.1)  the  problem  still  poses  a 
number  of  numerical  complexities,  primarily  because  of  ihe  required  limit  as 
.  However,  as  discussed  above,  (5.1)  does  provide  some  hope  in  solving 
what  is  otherwise  a  totally  untractable  problem.  It  is  tolt  there  is  suf¬ 
ficient  hope  as  to  warrant  further  investigation  of  the  utility  of  (5.1)  in 
the  solution  of  this  problem. 
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6.0  A  NEW  APPROACH  TO  COHERENT  SCATTERING  FROM  A  PERFECTLY  CONDUCTING 


RANDOMLY  ROUGH  SURFACE 
6.1  Background 

Among  those  Involved  with  the  applicetlons  of  rough  surface  scattering 
theory,  the  statement  Is  frequently  made  that  coherent  scattering  Is  reason¬ 
ably  well  understood  and  adequate  models  exist  for  the  phenomenon.  Because 
of  the  paucity  of  electromagnetic  scattering  data  [1],  one  must  go  to  the 
acoustic  field  to  appreciate  just  how  truly  erroneous  this  statement  is!  The 
acoustic  data  [2,3]  show  that  for  scattering  from  an  agitated  water  surface 
all  models  are  accurate  for  small  Rayleigh  roughness  parameter  .  However,  as 
either  the  frequency  or  surface  roughness  is  increased  or  the  angle  of  inci¬ 
dence  la  decreased,  the  data  show  a  significantly  stronger  scattered  field 
than  is  predicted  by  physical  optics  and  the  inclusion  of  shadowing  in  the 
model  only  makes  the  situation  worse  [4,5].  Tlie  acoustic  experiments  are  im¬ 
portant  because  they  were  designed  in  such  a  manner  as  to  eliminate  one  postu¬ 
lated  reason  for  why  early  electromagnetic  data  did  not  agree  with  the  physical 
optics  model  [6].  A  model  based  upon  pure  geometric  optics  has  been  developed 
[7]  but  it  tends  to  overestimate  the  mean  scattered  field.  Furthermore,  this 
analysis  appears  to  be  based  upon  a  questionable  transition  from  a  single 
sinusoid  surface  to  a  random  surface  and  it  provides  no  justification  for  the 
use  of  geometric  optics  for  a  situation  which  is  clearly  fraught  with  diffrac¬ 
tion  and  multiple  scattering  effects.  DeSanto  [8]  has  formally  solved  the 
problem  through  the  use  of  a  diagram  expansion  method  for  calculating  the 
stochastic  Green's  function  for  the  rough  surface.  DeSanto's  results  became 
even  more  significant  when  he  recently  showed  [9]  that  the  fii'st  correction 
term  to  the  r'lysical  optics  result  did  Indeed  increase  the  level  of  the  aver¬ 
age  scattered  field.  Unfortunately,  it  is  difficult  to  interpret  the  physical 
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ba9ls  of  the  higher  order  correction  terms,  each  of  whjich  Involve  the  solu¬ 
tion  of  an  Integral  equation  whose  complexity  Increases  with  order.  The 
need  for  further  analytical  and  experimental  research  on  this  problem  Is  there¬ 
fore  still  great. 

The  purpose  of  this  section  Is  to  present  a  new  approach  to  the  coherent 
scattering  problem  based  upon  averaging  the  magnetic  field  Integral  equation 
describing  the  current  Induced  on  a  perfectly  conducting  surface  by  an  inci¬ 
dent  field.  The  motivation  for  this  return  to  fundamentals  Is  as  follows. 

First,  It  Is  desirable  to  Investigate  solutions  to  stochastic  scattering  prob¬ 
lems  which  do  not  require  an  arbitrary  closure  assumption.  Second,  It  Is 
absolutely  essential  to  have  a  solution  wherein  mathematical  simplifications 
can  be  put  into  one  to  one  correspondence  with  physical  approximations.  Finally, 
acquiring  a  better  understanding  of  the  coherent  scattering  problem  Is  vital  to 
the  accurate  modeling  of  rough  surface  multipath  effects. 

6.2  Analysis 

The  analysis  will  be  presented  in  two  phases.  In  the  first  phase  the 
surface  roughness  will  be  assumed  to  be  arbitrarily  distributed.  In  the  sec¬ 
ond  phase,  the  surface  roughness  will  be  assumed  to  comprise  a  Gaussian  pro¬ 
cess.  Restriction  of  the  problem  to  a  Gaussian  surface  permits  the  detailed 
examination  of  certain  simplifying  assumptions  and  also  the  comparison  with 
De Santo's  [10]  results. 

6.2.1  Preliminaries 

The  rough  surface  Is  assumed  to  be  perfectly  conducting  and  Infinite  In 
extent.  The  surface  roughness  C(x,y)  Is  stipulated  to  comprise  a  zero  mean 
statistically  homogeneous  process  with  the  mean  surface  equal  to  the  z  0 
plane.  In  the  following  development,  position  vectors  will  be  denoted  by 

71 


t 


r  “  r^+zz  with  r^«xx+yy  and  for  a  point  on  the  surface  z*t(x,y). 

Using  an  exp(ju)t)  time  convention,  the  incident  magnetic  field  is  given  by 

H^(r)  -  hexp(  -  J  k^‘r) 

A  ^ 

where  for  vertical  polarization  h*y  ,  while  for  horizontal  polarization 
h  *  sin  9^  z  -  cos  O^x  and  the  karat  symbol  denotes  a  unit  vector.  The  field 
is  assumed  to  be  incident  along  the  positive  x-axis  so  the  incident  azimuth 
angle  is  also  zero,  "  0  •  The  Incident  wavevector  is  given  by 


k.  =  -k  (sin  0,  X  +  cos  9.  z) 
i  o  1  1 


where  9^  is  the  angle  measured  from  the  z-axis  or  the  normal  to  the  mean  sur¬ 
face  and  k  =  2Tr/X  is  the  free  space  wavenumber, 
o  o 

The  current  J  induced  on  the  surface  S  bv  the  incident  magnetic 

S  O' 

field  must  satisfy  the  magnetic  field  integral  equation  (MFIE) ,  i.e. 


J  (r)  = 


+  2^  fi(?)  X  xV^g(|^-"?^|)dS^ 


(6.1) 


for  t  •  In  (6.1)  n(r)  is  the  upward  directed  unit  normal  to  the  surface 
arjd  g(|r-r^|)  is  the  free  space  scalar  Green’s  function  where,  in  expanded 


form, 


^  -?:x-c9  +  z 

n(r)  =  - - 

yr  + 

’  X  y 


r  -  r  ) 
o 


exp(  -  j  k^|r-?^|) 


1^  ->• 

r  -  r_ 


It  should  be  noted  that  the  gradient  operating  on  g  in  (6.1),  when  evaluated 


on  the  surface  z  »  C  (x  ,y  )  ,  treats  the  random  height  as  if  it  were  inde- 
o  o  o  o 


pendent  of  the  coordinates  x  and  y 

o  o 


Expanding  the  double  cross  product 
in  (6.1),  converting  the  surface  integration  to  an  integration  over  the  z  *0 


plane  through  dS 


2  2  -> 

+  5  dr^  ,  and  multiplying  both  sides  of 


(6.1)  by  +  Cy  yields  the  following? 

J(r)  -  2n(r)  x  H^(r)  J’  |[n(r)*V^gl  J(r^)  -  [fi(r) •  J(r^) ]V^g  |  dr^  (6.2) 


where 


3<?)  - 


(6.3a) 


n(r)  «  -  -  CyY  +  2 


(6.3b) 


and  the  integration  is  over  the  entire  z^  «  0  plane.  For  future  reference, 
the  quantity  ?(r)  will  be  called  the  equivalent  flat  plane  current  because 
it  is  referenced  to  the  z  ■  0  plane.  Using  the  fact  that  ^  (r  )  must  be 

I  2  2 

tangential  to  the  surface  and  +  r,  >0,  there  results 


J(r)“C  J(r)+C  J(r) 

z  o  X  X  o  y  V  o 

o  •^o  ^ 


(6.4) 


Equation  (6.4)  can  be  substituted  in  the  right  side  of  (6.2)  to  yield  coupled 

Integral  equations  for  J  (r)  and  J  (r)  .  The  coupling  is  a  consequence  of 

X  y 

the  term  [n(r)»J(r  )]V  g  which,  with  the  substitution  of  (6.4)  in  (6.2) 

^  0  0 

yields  the  following  x  and  y-components; 


§ 

All  limits  on  the  Integrals  in  this  section  are  (-«>,<»)  so  they  will  not  be 
explicitly  shown. 


)J  (r  ) 
y  o 


} 


Jg. 


where  fi  is  x  or  y  .  It  should  be  noted  that  these  terms  are  proper- 
0  0-^0 

tionaltothe  difference  in  slopes  at  r^.  and  .  Thus,  if  the  surface  is 

o 

very  gently  undulating,  these  terms  should  be  very  small.  In  the  analysis  to 
follow,  these  terms  will  be  ignored;  titus,  the  problem  reduces  to  the  assump¬ 
tion  of  no  depolarization  or  the  case  of  a  surface  having  slopes  which  are  very 
slowly  varying  with  r^  .  Actually,  the  analysis  can  be  carried  through  for 
the  vector  or  coupled  equation  problem  in  essentially  the  same  manner  as  to 
be  presented  here.  However,  because  it  does  tend  to  symbolically  complicate 
the  equations  it  is  better  to  introduce  the  approach  with  the  scalar  problem. 


6.2.2  Arbitrarily  Distributed  Roughness 

Ignoring  the  n^J  term  in  (6.2)  yields 


j/r)  = 


2q*[fi(r) 


'*'1 

xH^r)] 


J-  f  f.r 

2iT  J  ^  S 


-Ifi- 


} 


J  (r  )  dr 
q  o  t 


(6.5) 

where  q  =  x  or  y  .  Computation  of  the  equivalent  flat  plane  current  is  not 
truly  Che  desired  end  result;  what  is  really  sought  is  the  average  scattered 
field  which  in  the  Fra'jnhofer  zone  is  proportional  to 


J  (r)  exp 


ak^^r.)>exp[j(k^^x  + 


(6.6) 


(where  it  has  been  assumed  that  the  averaging  operation  denoted  by  <•>  and 

the  surface  integration  can  be  interchanged).  The  averaging  operation  in  (6.6) 

implies  an  average  over  C  and  all  other  random  variables  upon  which  J(r) 

depends.  Clearly,  from  (6.5),  .r(r)  depends  on  the  slopes  C  and  C  5  further- 

X  y 

-y 

more,  experience  indicates  that  J(r)  should  also  depend  upon  the  curvature 


components  of  the  surface  (C  ,£;  ,C  )  •  In  point  of  fact,  J(r)  evaluated 

XX  xy  yy 

at  depends  upon  all  higher  order  derivatives  of  the  surface  heiglit  which 

are  correlated  with  the  surface  height  and  slopes  ev'tluatcd  at  r^  .  This 

o 

means  that,  for  the  general  case,  J  depends  upon  aU.  higher  order  derivatives 
of  the  surface  height.  Thus,  in  order  to  accomplish  the  average  in  (6.6)  the 
quantity  J(r)exp(Jk  ?)  must  be  multiplied  by  the  single  point  .loint  pro)^- 

8Z  ^ 


billtv  density  function 


Pj(C,V^,7^r,,y^C,-..) 


2 

and  averaged  over  C»VC,7  ?,•••»  where  7  4  is  a  symbolic  notation  for  all 


th 


order  derivatives  of  the  surface  height  evaluated  at  the  point  r^  .  That 


is. 


<J(r)exp(jk  0> 

9  Z 


SH 


J(r)exp(jk^j^  C)  pj^(C,Vi;.V 


•  dCdV?*-* 


(6.7) 


The  right  side  of  (6.7)  can  also  be  written  as  the  convolutions  of  the  infinite 
dimensional  Fourier  transforms  of  J(r)  and  p^(*)  as  follows; 


<•>  ■  lim 

n-xxj  (2^) 


•  dBj^  d&2  ^33*  • 


(6.8) 


where  B  is  an  n-dlmensional  'Vector"  and 
n 


•  dC  d7c---d7"”^C 


(6.9) 
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(6.10) 


‘tj^Ckgz  "&l«"^2’' ■ ‘’"®n^  “ 


- J  (32  •^^+-  •  •  -^Sn •^''’^^)]  dVC*  •  •  dV’'"^C 


The  notation  Is  a  hit  cumbersome  here,  co  explicitly  writing  a  few  terms  may 
be  helpful,  e.g. 


@2 -VC  -  S,?^  +  By?, 


6,-V  ;.B?  +BC  +BC 

3  XX  XX  xy  xy  yy  yy 


6,  *7  c  =  3^  ■•■3  c  ■•■3C  ■••3C 

4  XXX  XXX  xxy  xxy  xyy  xyy  yyy  yyy 


Since  <I)^  is  the  Fourier  transform  of  Pj^,  it  is  the  single  point  loint  charac- 

<K/ 

terls  ’  fxmction  for  the  random  surface.  Note  also  that  J  is  the  Fourier 
transform  of  J  with  respect  to  all  the  random  variables  upon  which  it  depends 
(an  iaflnite  number  in  general). 

According  to  (6.8),  J  is  required  in  order  to  compute  the  average  scat¬ 
tered  field.  This  suggests  (6.5)  should  be  multiplied  by 

OO 

exp(J  ■f  j  ^ 
n=l 


and  then  averaged.  By  expressing  the  averaging  integrations  as  convolutions 
of  Fourj.er  transforms,  an  integral  equation  for  J  in  (kj^,k2»k2»*  *  *  )-space 
can  be  obtained.  The  average  of  the  term  on  the  left  side  of  (6.5)  is  given 
by 
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■  h';. . 


<J^rt)exp(jkjC+J  ^  2  k^^j 

n»l  *'  n=l 


«v"o 


•  Pja.VC,-*- )dc  dV^- 


llm 
n-H»  (27r) 


de^  de2 


(6.12) 


The  source  term  cn  the  right  side  of  (6.5)  can  be  written  as  follows; 


2^-[nXH^(^)]  =  2H^[C^+C^^^  +  CyCy]  exp  (-jk^^C- jk^4^)  (6.13) 

where  C  ,  C  and  C  are  determined  by  the  polarization  of  H  but  are  inde- 
o  X  y 

pendent  of  the  random  surface  variables.  Thus,  the  average  of  tne  product  of 
this  term  and  (6.11)  is  given  by 

x;i(J))>  .  2H„  exp  H'  ] 


2y 


<l>l(ki-ki^,k2,k3,*--) 


(6.14) 


where  k2  is  symbolic  for  the  variables  k^^  and  k^^  .  The  (27r) 

term  does  not  appear  in  (6.14)  because  no  convolutions  are  required  -  only 

straightforward  Fourier  transforms.  The  symbols  and  3.  denote  the 

2x  2y 

partial  derivative  operators  3/3k2^  and  respectively. 

The  average  of  the  product  of  the  Integral  term  on  the  right  side  of  (6.5) 
and  the  exponential  factor  in  (6.11)  is  somewhat  more  involved  than  the  aver¬ 
ages  of  the  other  terms  in  (6.5).  First,  the  two  point  Joint  probability  den¬ 
sity  function 


mmmmmmmmrnmmmmmmmim 


.1  It kukAk  i . 
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2  2 
w  r-»^ 


"V  • 


Po(C»C  .VC.V  c  ,V  C.V  5  ,  •••  ;  Ar.  •  -  r.  ) 

2o  oo  oo  tct 


must  be  used  because  of  the  additional  set  of  random  variables  C 

o  o  Of  o  o 

in  this  term.  Thus,  the  average  can  be  written  as  follows; 

(  J(-  ^  ^  ^  1 


•  P2<C,5^,V;.V^r,^,  •••  )  dr^.  dr.  dC^  dVC  ••• 

o 


(6.15) 


Assuming  that  the  order  of  the  integrations  can  be  arbitrarily  interchanged, 
the  C~lntegratlon  can  be  written  as  a  convolution  of  the  ^-Fourier  transforms 
of  the  Green's  function  derivatives  and  P2(*)  •  Noting  that 


Fj  I - 9^^ - 1  -  -  exp  F^ 


9g(rj.-r^  ,0 


where  denotes  the  Fourier  transform  with  respect  C  and  is 


transform  variable,  and  substituting 


g(Arj.,3^)  =  F^  |g(Ar^  ,?;)  I  =  J  g(Ar^,5)exp(j3^C)di; 
^  f  3g(Ar  ,C)^  r  9g(Ar  .C) 

l^(Ar^,3o)  =F^  [ - j  =J  - -  exp(j3^ 


9g(Ar  ,0 

- rr -  exp(j3„C)dC 

O 


in  the  convolution  integration  with  Ar^  =  r^ 
for  (6.15); 


-  r^  yields  the  following  form 


3g(Ar  ,3^)  3g(Ar  ,3  ) 

r  to  ^  to  ^  /A^ 

3x  “  S  3y  ' 
o  ^  o 


f8o>) 


CO 

2  ®n+l‘''”')P2'V®o-^o’’'’''<.5o’"‘>‘‘Pt  '*®o‘‘ V''MVo' ■  ■ 


(6.16) 


where  the  tilde  symbol  denotes  the  Fourier  transform  of  P2(*)  with  respect 


to  ? 


The  Integration  in  (6.16)  can  also  be  written  as  a  convolution  as  fol¬ 


lows; 


3g(Al  ,3  )  3g(Ar  ,i:  )  ~  h- 

axV  -  S  . - 


00 

•  exp  (j  ^ 

^  n-1  ^ 


o 

(6.17) 


where  the  tilde  over  J  and  the  second  tilde  over  p_  denote  the  Fourier 

q  ^ 

2 

transform  with  respect  to  C  •  The  remaining  Integrations  over  V^.V  C.*'* 

o 

are  simply  Fourier  transforms,  so 


3g(Ar  ,3„)  3g(Ar  ,3^)  ^ 

- \  J - 3^  \  ' 

"^o  4x  ^o  ‘^2y  ^ 


\^^o»®l^  P2^V^o’V®l’^2’^o^o’S’^o^o’  ’ 


••)  dl  d3  d3,dV^C  dvjc  '** 

t  o  1  o  o  o  o 
o 

(6.18) 


where 


p  (k  -0  .3-0, ,k„,V  c  ,k  • *) 

^2  1ool2oo3oo 


•  dVcdV^C*’ 


2.  r,3. 


is  the  Fourier  transform  of  p„  over  all  variables  except  7  C  ,7  C  .7  C  ,*•* 

2  - ^  o  o  o  0  0  o 

It  should  be  noted  that  the  differential  operators  3,  and  3,  in  (6.18) 

‘"2y 

operate  only  on  p^  •  The  integrations  in  (6.18)  over  ^o*’o*'^o  ^o’^o  ^^o’ *  ' 
can  be  written  as  convolutions  of  the  Fourier  transforms  of  J  and  p*  ,  i.e. 

q 


<  •  >  =  1  im 


n-x»  (2Tr)(27i) 


1 _  f  ff  sKA^.e,) 

(2,)"'  J'JV  K 


3g(Ar  ,B„) 

\  *  J  — \ 

^2x  ^o  *^2y 


3  (r^  ♦3j^.32»^3****)P2  ’^o~^1  *^2 ’"^2 *^3 ’”^3* *  ’  * *^^1* 

o  «  o 


(6.19) 


where 


o 


n=l 


•  dV  r,  dV  c  ’  •  • 
00  00 


P-,(k  -B  ,3  -3,,k-,(?„,^~,^»,'"')  =  I  p,(k  -3  ,3  -3  ,k_,7  ?  ,7^5 

-2  10012233  J  J^2  1ool2oo3oo 


00 

(j  2  "V-o 


dV  c  ’** 

o  o 
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The  function  Is  the  transform  of  the  two  point  .joint  probability  den- 

sity  function;  thus,  it  is  equal  to  the  two  point  joint  characteristic  func¬ 
tion  <t>2  •  Substituting  (6.19).  (6.14),  and  (6.12)  back  into  the  average  of 
the  weighted  [by  the  term  in  (6.11)]  equation  (6.5)  yields 


2H  exp  (-jl  ‘r  j[c  +  j  C  3  +  j  C  3  1  (J) 

0  1  t  [_  o  X  k2^  y  K2yJ 


1  ^ 


+  lira 


n-*«  (2Tr)  (27r) 


1  r  r  r .  ,  is. ,  . ; 

^J-J  -Tj;- ^  S, 


(6.20) 


Substituting  Ax  =  x  -  ,  Ay  =  y  -  y^  ,  and  Ar^  =  ”  r^.  in  the  right  most 

term  in  (6.20)  yields  the  following; 


n-»«  (2tt)‘ 


\  ^'^t  •  ^1  ’^2  ’  *  ’  *  ^  '*^1  ’^2"^2  ’  ■  ■ '  ^  '^^1  *^^2 


exp  ^  ^x\,,  ^  ^  ^^^2  ^ 


+  lim 
n-x» 


I  3i(Ar  ,3  )  ^  ^ 

-3,  +  j  - aX -  K  ■*■  S^CAr  ,3  )> 

\x  2y  C  t  oj 


‘*'2^V®o’V^1’^2’"^2’^3’“^3’'“^  ?q  ^ V^^t ’^1  *^2 ’^3’ ’  *  *  ‘*^o*^®l‘*^2‘ 


The  single  poLnL  Joint  cliaracterlstic  function  is  Independent,  of  Ar^ 

wliile  the  two  point  Joint  characteristic,  function  depends  on  through 

the  correlation  function.  Since  all  terms  in  (6.21)  must  exhibit  the  same 
dependence  upon  and  only  J  and  exp(-jk^*r^)  are  functions  of  r^  , 


(6.21)  implies  that  can  be  written  as  follows; 


'^q^^t’^1’^2’*  *  ’  ^  =  jq(3j^.32.“  *)  exp  (-jk^*r^) 


(6.22) 


Substituting  this  result  in  (6.21)  and  rearranging  terms  produces  the  following 

->■ 

integral  equation  for  Jq^3j^.32* '  ‘  *  )  ♦ 


11m - - 

(2Tr)' 


~„T  [j-  •  •/  >  [*1  ('‘r<*cM2  • 


■  (2.7  J  J  V  '  3^-  ■  ■  ‘‘2X  ^  '‘2y  ] 

■  2“o['=a  3  J  ‘V‘ix4-*3'---> 

Before  a  detailed  discussion  of  (6.23)  is  presented,  it  is  advisable  to  review 
the  mathematical  meaning  of  the  various  terms,  e.g. 

.iq(3i.32.  •  •  •  )  =  exp(jit^*r^)  'J"  d^(r|.,f,,VC,V^?;,*  •  •  )exp  ^  ^  2  ^n+1^"^  ] 


•  dCdVc  dV  C  • 


(6.23a) 


'/"A 


)expj^. 


J(kj-e^)C-hJ 


n*=l 


•> 

n+l  ^ii+l  ^ 


/"c] 


dCdVr;  dV  c*“ 


g(A^t»Bo)  =  J' g(Arj.,r,)  exp  (J3^C)d^ 
r  3g(Ar  ,C) 

t’^o^  ‘  J  - 3C - exp(j3/.)dC 


g^CAr 


(6.23b) 


(6.23c) 


(6.23d) 


and 


e*p[jOCi-6„)C  +  J  f  -  i  Y 

n=l  n=l 


•  dC  dC^  dVC  dV^^^  •  *  •  (6.23e) 

Also,  equation  (6.22)  confirms  DeSanto's  earlier  analysis  [8]  in  that  it  shows 
that  the  coherent  scattered  field  is  specular  in  nature.  That  is,  if  the  aver¬ 
age  in  (6.6)  is  written  as  a  convolution  such  as  in  (6.8)  and  (6.22)  is  sub¬ 
stituted  for  ,  the  r^-spatial  integration  will  yield  a  product  of  6-func- 
tions,  e.g. 


/■ 


(6.2A) 


which  shows  that  the  scattered  field  is  nonzero  only  for  k  =  k.  and 

sx  ix 
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6.2.2. 1  Discussion  of  Results 


Equation  (6.23)  is  an  integral  equation  for  this 

function  can  be  determined,  the  exact  amplitude  of  the  scalar  scattered  field 
may  be  computed  as  follows; 

.  •  •  • )  ■  ■  >“»!  ■'^2  '‘®3'  •  • 

(6.25) 

The  quantity  from  (6.22),  is  proportional  to  the  Fourier  transform  of  the 

equivalent  flat  plane  current  with  respect  to  all  the  random  variables  of  which 

this  current  is  a  function.  Caution  should  be  exercised  in  any  attempt  to  at- 

2 

tach  a  physical  meaning  to  because  in  transforming  from  C»'** 

to  ,32  stochastic  character  of  the  random  variables  is  lost. 

In  fact  this  is  the  fundamental  reason  behind  working  in  3'j^,32  »^3»  *  ’  *  space. 
Hiat  is,  if  one  were  to  average  (6.5),  weighted  by  the  exponential  factor  in 
(6.11),  in  the  conventional  manner  of  multiplying  by  the  appropriate  joint 
density  functions  and  integrating  over  all  random  variables  directly  without 
going  to  the  transform  space,  the  stochastic  nature  of  the  random  variables 
would  prohibit  one  from  obtaining  a  single  equation  such  as  (6.23).  It  is 
well  known  that  a  conventional  average  of  integral  equations  such  as  (6.5) 
leads  to  an  infinite  s^  of  integral  equations  [11]  because  one  does  not  know 
the  average  of  the  pradu.ct  of  and  the  kernel  inside  the  integral  in  (6.5). 

VAjat  lias  been  shown  here  is  that  if  the  averages  are  expressed  as  convolutions 
in  l  rails  form  space  rather  than  direct  integrations  over  the  random  variables, 
it  is  possible  to  obtain  a  single  integral  equation  because  the  transformed 
product  of  .1^  and  the  kernel  term  can  be  factored.  The  price  that  one  pays 
for  the  single  integral  equation  is  that  it  has  infinite  dimensions  because 
all  order  derivatives  of  which  have  a  nonzero  two-point  correlation  with  5 


lim - - — r  I 

(2n)*^  J 
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must  be  Included.  In  fact, the  primary  difference  between  this  approach  and 
the  conventional  method  [11]  which  gives  rise  to  an  infinite  set  of  equations 
Is  the  following.  In  the  conventional  approacli,  one  attempts  to  solve  for 
the  average  of  the  desired  unknown  quantity  without  trying  to  explicitly  deter¬ 
mine  its  dependence  upoti  the  randcn  parameters  in  the  problem.  In  the  approach 
presented  here  the  exact  opposite  is  done  only  in  tlie  transform  domain  where 
the  stochastic  character  of  the  random  parameters  enters  only  through  the 
kernel  and  the  correlations  between  parameters.  Furthermore,  in  the  conven¬ 
tional  approach  it  is  frequently  difficult  to  attacli  physical  signific.ance  to 
auy  truncation  (closure)  or  partial  summation  of  the  infinite  set  of  equations. 
In  the  method  presented  here,  truncation  of  the  infinite  dimensionality  of 
(6.23)  is  determined  entirely  by  the  relative  magnitude  of  the  correlations 
between  the  random  variables. 

There  is  one  notentlai  problem  with  (6.23)  which  may  make  it  less  attrac¬ 
tive  than  the  conventionaJ.  approach.  In  the  Infinite  equation  solution,  all 
the  integral  equations  are  of  the  second  kind  and,  thus,  normally  amenable 
[10]  to  numerical  solution.  On  the  other  hand,  (6.23)  is  of  the  first  kind 
which  is  usually  rift  with  problems  [12].  This  is  certainly  a  point  to  be 
considered  in  the  future;  however,  there  are  two  reasons  why  it  may  not  be  a 
problem.  First,  the  desired  quantity  is  (6.25)  and  not  simply  ;  thus 
problems  associated  with  the  accurate  recovery  of  from  (6.23)  may  vanish 

when  computing  (6.25).  Second,  since  (6.23)  is  obtained  from  an  integral 
eqviatlon  of  the  second  kind,  this  may  also  minimize  some  of  the  problems  nor¬ 
mally  associated  with  equations  of  the  first  kind. 

There  are  a  number  of  interesting  results  that  can  be  obtained  from  (6.23) 
without  specifying  the  forms  of  the  one  and  two-point  joint  characteristic 
functions.  In  the  physical  optics  approximation,  the  term  in  (6.23)  involving 


H5 


■Iml/  li.i.  .  ......  <  It  . . . 


the  Green's  function  is  ignored  and,  according  to  (6.25),  the  aicplitude  of  the 
average  scattered  field  is  thus  given  by 


k„  -H) 

2y 


The  terms  involving  the  partial  derivatives  with  respect  to  k2^  and  k^ 


are  equivalent  to  the  following  average; 


(6.26b) 


where  q  =  x  or  y  .  For  a  statistically  homogeneous  process,  C  is  uncor¬ 
related  with  3C/3x  and  Sf^/Sy  [1)1.  Consequently,  the  average  in  (6.26b) 
is  the  product  of  averages  of  the  slope  and  height  factors  and  if  the  slopes 
are  a  zero  mean  process  then  (6.26b)  is  identically  zero.  Thus,  (6.26a)  re¬ 
duces  to  2H  C  (J)  (k  -k.  )  where  <[>(♦)  is  the  marginal  characteristic  func- 
tion  for  the  surface  Ixeight.  In  the  physical  optics  limit  the  specularly  scat¬ 
tered  field  is  independent  of  the  slopes  provided  the  surface  is  statistically 
homogeneous  and  the  slopes  are  zero  mean. 


The  function  corresponding  to  the  physical  optics  approximation 

(denoted  by  )  may  be  obtained  from  (6.23)  by  inspection,  e.g. 


-  .1C^«(p2^)6'(S2y)J  (6.27) 

here  6(*)  is  the  Dirac  delta  and  <S'(»)  is  its  derivative.  The  next  level 
of  approximation  is  the  so-called  Bom  approximation.  In  this  approach,  one 
solves  for  ~  ^^s^  moving  the  term  containing  the  Green's  function 
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trsnsfonn  to  the  right  side  of  (6«23)  and  substituting  [from  (6.27)1  In 

this  term.  This  leads  to  an  approximation  for  the  complex  amplitude  of  the 
scattered  field  which  Is  valid  wlien  the  term  containing  the  Green's  function 
transform  is  small  compared  to  the  source  term;  In  particular,  the  result  is 
as  follows 

<e/  .  2H„  llm  [C^  +  J  0^8  +  1 


2H  Ilf 

■*  2  I  I  §Ax 

(2^)  a^x'Sx-^^' 


r  fr  3g(Ar  .6  )  3g(Ar  ,B  )  ^  ^ 

iL°- — 5i— 


e2y‘‘2y-^“ 


■5  "J  ^ 

2  sz  o  o  Iz 


X  33 


2x 


-  J  c 


^2  (*^sz-^ 


y  33 


2y 


,3  -k.  ,0,k-  .0,-3»^,0,0,--*)'\  ^ 

o  0  iz  2y  2y  /  dAr  d3  (6.28) 

J  to 


where,  as  a  reminder,  the  arguments  of  <t>2  at-e  the  transform  variables  cor¬ 
responding  to  the  following  order  of  random  variables  (C»Cjj*3c/9x,3C/3y , 

3C  /3x  ,3c  /3y  .V^C.V^C  ,,*••) .  The  B  superscript  on  <E  >  in  (6.28)  denotes 

*  O  *  ^  O  O  c> 

the  Bom  approximation. 

One  final  result  that  should  be  demonstrated  is  the  limit  of  a  perfectly 
flat  plane.  In  this  limit  the  correlation  function  for  the  heights  goes  to 
the  mean  square  height,  <C  =  a  «  <C  >  ,  while  the  mean  square  slope, 

curvature,  rate  of  change  of  curvature,  etc.,  go  to  zero.  Consequently, 


..i. 


f 


'1 
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'Che  Green’s  function  transform  g,(Ar^,3  )  Is  an  odd  function  of  3  because 

G  t  o  o 

g>  “  J  ®nd  g  is  an  even  function  of  3  •  Thus,  the  term  in  (6.23) 

>  o  o 

involving  the  Green’s  function  transforms  Is  zero  and  (6.23)  reduces  to 

which  yields  the  proper  value  for  <E  >  as  k,->k  and  also 

S  X  s  z 

j^(3,  =  2H  C  6(3,-k.  )  n  6(3.)  (6.30) 

q  1  2  J  o  o  1  iz  £_2  i 

Furthermore,  the  inverse  transformer  (6.22)  is  given,  in  this  case,  by 

Jq(r^,r,)  =  2H^C^exp(-j  k^^q-jt^*rj.)  (6.31) 

and  tills  is  a  valid  transformation  since  C  does  not  depend  upon  r^  .  It 
siiould  be  noted  that  tliis  is  the  case  of  a  randomly  elevated  plane  and  yet  it 
produce.s  the  same  results  as  the  physical  optics  approximation.  Thi'^  corre¬ 
spondence  leads  to  increased  suspicion  of  the  physical  optics  result.  That 
is,  physical  optics  treats  tlie  problem  as  a  randomly  elevated  plane. 

6.2.3  Jointly  Gaussian  Distributed  Roughness 

The  results  of  the  previous  section  are  important  because  they  provide  a 
rigorous  mathematical  foundation  for  the  multivariate  approach.  They  are  also 
general  in  that  they  are  valid  for  any  zero  mean  statistically  homogeneous 
rougliness.  Wliile  it  is  possible  Co  obtain  certain  asymptotic  solutions  such 
as  witli  the  physical  optics  and  Born  approximations  or  In  the  case  of  a  ran¬ 
domly  elevated  plane,  it  is  difficult  to  appreciate  the  power  of  this  approach 

without  seeing  it  applied  to  a  specific  surface  height  distribution.  This  is 
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lim 
n->oo  (2iT) 


particularly  true  In  regard  tu  reducing  the  dimensionality  of  the  integral 
equation.  That  is,  given  the  form  of  the  single  and  two  point  joint  char¬ 
acteristic  functions  and  (6.23),  how  does  one  go  about  determining  which  sur- 

2 

face  parameters,  l.e.  C«VC,V  important  and  over  what  range  of 

values?  To  accomplish  this  goal,  the  jointly  Gaussian  surface  has  been  se¬ 
lected.  There  are  two  reasons  for  this  choice.  First,  the  Gaussian  surface 
has  been  extensively  studied  by  others  (2-101.  Second,  the  jointly  Gaussian 
density  and  characteristic  function  have  known  closed  mathematical  forms; 
something  which  is  difficult,  at  best,  to  obtain  for  other  distributions. 

The  purposes  of  this  section  are  to  demonstrate  how  one  goes  about  solving  (6.23) 
for  a  Gaussian  surface,  obtain  asymptotic  solutions,  ai^d  compare  these  results 
with  others. 

The  surface  is  assumed  to  be  zero  mean,  jointly  Gaussian,  and  statis¬ 
tically  homogeneous.  With  U2  the  column  matrix 


C 

^o 

VC 

V  c 

o  o 


(6.32a) 


and  C2  the  square  covariance  matrix 


<C^> 

<C’Vc>  •• 

A 

0 

V 

<VC*C> 

<(VC)^>  • 

<v  c  •c> 

00 

A 

< 

0 

0 

• 

0 

V 

<7„c„-Vr.> 

(6.32b) 
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^,»^;/rt^«.  -.  •^-•..'i^v..  :^  ,,rs^  I r>  'fiw tt w>.‘«r*- W -  -* ».Wi 


-T  ,=  x 

h  <'=2> 


"2A] 


the  two  point  joint  probability  density  function  is  given  by 

=  lira  [(^'")”  1^2!]  exp  ["  ^2^  (^2) 

(6.33) 
—  T 

where  the  T  superscript  denotes  the  transpose  of  U2  t  i*e.  U2  is  a  row 
matrix,  and  1^2!  *^he  determinant  of  the  covariance  matrix.  With  V2 

the  column  matrix  of  transform  variables 


V  -  k 

2  ^^2 


(6.3Aa) 


the  two  point  joint  characteristic  function  is  given  by 


<6. 34b) 


The  terms  in  the  covariance  matrix  are  typically  as  follows i 


<C  l,>  =  R(Ar^) 
<r  r,^>  =  0 

<r,  r,  >  =  0 

y 


-R  (ArJ 

yo 

y  t 
->• 

<K 

c  >  - 

R 

(ArJ 

XXo 

XX 

t 

<c 

C  >  = 

R 

(ArJ 

yyo 

yy 

t 

H 

A 

VJ 

-R  (ArJ 

X 

XX  t 

<C 

C  >  = 

-R 

y 

yo 

yy  t 

•4- 

<c 

X 

c  >  = 

XXo 

R  (Ar.) 
xxx'  t 

<C 

C  >  = 

R  ,„(Ar.) 

X 

yyo 

xyy  t 

c  >  = 

R  (Ar  J 

y 

XXq 

yxx  t 

<'Q 

C  >  = 

R  ,(Ar^) 

y 

yyo 

yyy  t 

where  Ar. 


=  r^  -  r,  ,  the  x  and  y-subscripts  denote  differentiation  with 


respect  to  Ax  and  Ay  ,  respectively,  and  K'Ar^)  is  the  surface  height 
correlation  function.  The  above  results  follow  from  the  general  relationship 

ri3] 

.m+n+p+q 


,,in-hi  9 

J — -k.  . _ 

.3x™3y"  3xP3y^^ 


(-1) 


p+q 


3Ax"^3Av"^ 


(6.35) 


where  m,  n,  p,  and  q  are  nonnegative  integers. 
If  Uj^  is  the  column  matrix 


and  is  the  covariance  matrix 


(6.36) 


1 

2 

> 

<r;*V^C>  *" 

c,  * 

<VC*r,> 

<('?C)^> 

<VC‘V^C> 

1 

<v^t;»vc> 

2  2 

<(S  O  >  "* 

i 

* 

- 

(6.37) 


the  single  point  joint  prebability  density  function  is  given  by 

P2(C,Vcy;,-'*)  =  lira  [<2Tr)"  |c^l]^^  exp  [-u^(:c^)  (6. 


n-H» 


and  the  corresponding  joint  characteristic  function  is 


where  is  tike  column  matrix  of  tran.siorm  varlahles 
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38) 


(6.39) 


. . . . .  » 

....... 


The  diagonal  elements  In  (6.37)  are  the  mean  square  (or  variance)  parameters 

2  2 

of  the  surface.  The  mean  square  height  Is  O  =  >  ,  the  x  and  y-components 


of  the  mean  square  slope  are,  respectively. 


X 


<C^>  , 
y 


one  can  similarly  define  the  mean  square  curvature  components,  rate  of  change 


of  curvature,  etc.  It  should  be  noted  that  these  variance  parameters  can  be 
obtained  from  (6.35)  evaluated  at  Ar^  =  0  .  Furthermore,  in  order  for  these 
variances  to  be  finite,  there  are  certain  analytic  properties  required  of 
R(Ar^)  in  the  neighborhood  of  Ar^  =  0  . 

Equations  (6.34)  and  (6.39)  provide  the  joint  characteristic  functions 
required  in  (6.23).  However,  these  functions  still  depend  upon  an  infinite 
number  of  transform  variables  and  this  will  make  the  meinipulations  tiaore  dif¬ 
ficult  CO  follow.  The  crucial  points  to  be  made  can  be  accomplished  just  as 
completely  if  only  two  transform  variables  are  considered.  Since  the  two  most 
important  surface  parameters  are  usually  the  height  and  slopes,  these  will  be 
selected  for  the  demonstration.  In  order  to  eliminate  the  other  surface 
parameters,  it  will  be  necessary  to  assume  a  very  gently  undulating  surface 

which  has  a  vanishingly  small  mean  square  curvature,  rate  of  change  of  curva- 

n  2 

ture,  etc.;  that  is,  <(V  C)  >  -  0  for  n  >  1  .  Under  these  conditions. 


45^  (k^-p.j  ,k2~P2  ’^3"^3’  ■  "  ^  ^  ^1  ^‘T^l  ’^2"^2^ 

%  ''‘^l“^'o*^o”^''l  ’'^2  ’”^2  '^2  ’®o~^l  ’*^2  ’"'^2^ 


(6.41a) 

(6.41b) 
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and  In  (6,23) 


r 


(6.41c) 


because  the  joint  probability  density  functions  go  to  Dirac  deltas,  6(V"c)  , 
while  is  Independent  of  the  •  Under  these  circumstances,  (6.23) 

reduces  to  the  following  form; 

I, >[*(4-6, //{>  4^^ 


(2tt)' 


+  j 


38(Ar  ,B  ) 


■t 


2“ol"o  J  i 


} 


(6.42) 


where,  because  the  height  and  slopes  are  independent  for  a  Gaussian  surface, 
the  single  point  Joint  characteristic  function  b  is  been  factored  into  a  prod¬ 
uct  of  marginal  cncracteristic  functions  for  the  height  and  slopes.  It  should 
be  noted  that  (6. <>2)  la  exactly  the  equation  that  would  result  from  assuming 
that  the  c'.;rfen>  '■"'ipenda  only  upon  the  random  surface  height  and  slopes. 

However,  as  noted  above,  (6.42)  is  valid  only  when  the  variances  or  mean 
square  values  of  all  the  higher  order  surface  derivatives  are  vanishingly  umall. 

The  remaining  portion  of  this  section  shall  be  devoted  to  two  goals. 

Using  the  height  and  slopes  as  examples,  an  attempt  will  be  made  to  obtain 
some  quantitative  measure  of  when,  say,  the  slopes  can  be  neglected  relative 
to  the  height.  Hopefully,  it  will  then  bo  possible  to  generalize  this  result 
to  the  point  of  also  estimating  the  required  smallness  of  the  curvature 
variance,  the  rate  of  cha  go  of  curvature  varlfince,  etc.  in  order  tiial  tliey 
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p',.  .  •  j  ■  ■  .  ... 


may  also  be  Ignored.  This  Is  an  extremely  imporrant  issue  alree  the  solution 
of  (6.23)  in  any  more  than  three  dimensions  (height  and  two  components  of 
slope)  is  highly  impractical.  The  second  goal  is  to  study  the  possibility 
of  solving  (6.42).  That  is,  since  (6.42)  represents  surfaces  which  are  not 
unreasonable,  its  solution  would  represent  a  considerable  breakthrough  in  the 
understanding  of  scattering  from  rough  surfaces. 

After  some  straightforward  manipulation  of  the  characteristic  functions, 
(6.42)  may  be  written  as  follows; 


///iq<erB2x-«2y>'*’<'‘r6l>*0'2x-e2x>4‘<V' 


i2vy 


R)(3  -3.)  (3  -k.)  -  (S  +R  )B,  k- 
o  1  o  1  X  XX  2x  2x 


(S  +R  )3,  k-  -  (k,-3  )(R  3.,  +R  3.,  )  +  (3  -3,)(R  k.  +R  k-  ) 
y  yy  2y  2y  1  o"^'  ^  2x  y  2y  o  1  x  2x  y  2y 


2y 


=  2H 


(6.43) 


2  2  2  2 

where  S  =  <?;>,  S  =<C>,  and  G  is  the  term  involving  the  Green’s 
X  V  y  y 

fund  Von  transforms  which  is  as  follows; 

->  9g(Ar  ,3  )  r  ,  3  38(^r,.»3„)  r-  , 

G(Ar^,3„,k2,(>.2)  =.i  -'-3aV-[-^^‘^2x-P2x>  “W"- [‘ V  V^2y> 


+  3. 


2y 


(6.44) 


The  term  in  (6.43)  involving  G  does  not  converge  very  rapidly  as  Ar^  » 
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so  to  avoid  potential  problems  this  asymptotic  behavior  will  be  subtracted 
out.  That  Is,  since  the  correlation  function  and  its  derivatives  go  to 
zero  as  Ar^  ,  this  terra  asyraptotically  approaches,  for  Ar^  sufficiently 


large. 


exp  [-  -s>2x'‘2x  -^%yW] 


exp(jic^‘Ar^)  dAr^  d3^  »  r(3j^,k^,’k^ 


(6.45) 


The  Ar^-lntegration  in  (6.43)  can  be  easily  accomplished  with  the  following 


result; 


•  (W  -®x»2x'‘2x  -  ^'®2y''2y] 

where  the  triple  tildes  denote  the  three  dimensional  (Ax,Ay,C)-l'''ourier 
transform,  the  subscripts  on  the  g  denote  the  derivatives  of  g(Ar^,C) 
with  respect  to  the  Indicated  variable,  and  k^^  =  k^^  x  +  k^^  y  is  the  trans¬ 
verse  part  of  the  Incident  wave  vector.  The  3^-lntegration  is  difficult  to 
perform  directly  because  of  the  Gaussian  factor  in  (6.46);  however,  it  can  be 
accomplished  by  using  Parseval's  theorem  to  rewrite  (6.46)  as  an  integration 
over  the.  product  of  inverse  Fourier  transforms  (with  respect  to  3^).  That  is, 
with* 


*The  4iT  factor  arises  because  g  is  4'n’  times  the  conventional  free  space 
Croi'ii's  ^unction  oxp(-.|  k^r)/4Tir  , 


'  ix 

-2T 


J.  ATTk.  expC-jk^^lBl) 

|(ku.,3^)exp(-j3^2)de^  "  ~2  Vr 


(6.47a 


~  -*■ 

- 2^ 


4TTk^^  exp(-jk^^lzl) 

)exp(-jB^z)d3^  - 


"  2Tr 


J 


3^)exp(-j3^z)de^  -  Y 


(6.47b 


(6.47c 


Kers  k  is  the  z-component  oi  the  incident  wave  vector,  and  [14,  pg.  63) 
iz 


1  z  >0 


sgn(z)  *  0  z  ” 0 

-1  z  <  0 


(6.47c 


A  J  exp  [-o^B„-8^)<3„-k^)  -  J  e^xJdS^  .  Q(x) 


(6.47 


(6.46)  may  be  written  as  follows; 


r(-)  =  2TTexp  [-^^^2xh^'^yW^2y 


1  /  ( 


S.^k.k, 


+  2it  sgn(z)  )  q(-z)exp(-jkj^lzl )dz 


(6.4£ 


Completing  the  z-integtat  Ion  yields 


r(-)  =  exp  [-  s\2xP2x  '  ^Sy®2y ] 


.iLArLui:..Li^^‘<.)l^v!aat’kx4iip. . i-.xui 
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where 


2  2 

+  exp  -o 


+  2,2e,p(-a\Bi)Lp[-a\2^[l-(!i^A]]]  .ertc[-^  (' ■(^■)} 


2  2 

exp  -  o 


(6.50) 


where  erfc(*)  is  the  complementary  error  fun. t ion.  Substituting  (6.49)  in 
(6. A3)  in  such  a  manner  as  to  regulariz.^  the  dAr^  integral  results  in  the 
following; 


exp(jkj^*Ar^)dAr^  d3^  -  T  M^l 


-  J  "x^x  W2  -  S/2J  <^-51) 


where  K(*)  is  given  below  and  the  addendum  of  "S”R*0,n“0,l,2"  means  that 
the  correlation  function  and  all  of  its  derivatives  are  to  be  set  t  ^  zero; 
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K(Ar,  “  (:(A,?^.f^^^,k.^.(y  exp  [-  (O^-K)  (P^-3^)  (3^-k^  > 


-(S^  +  R  )3»  k_  -(S,^+R  )  3,  k.  ~(k,  -3„)(R3,  +R  3o  ) 

'  X  XX  2x  2x  y  yy‘  2y  2y  1  o  '  x  2x  y  3y 


+  (3  -3,)(R  k^  +R  k,.  )  -  R  (3,  k,  +3,  k.  ) 
'  o  ^1''  X  2x  y  2y  xy  2y  2x  2x  2y 


] 


(6.52) 


Equation  (6.51)  Is  an  exact  Integral  equation  for  the  function  under 

the  conditions  that  all  surfree  height  derivatives  of  order  two  or  larger  are 
negligible  In  a  mean  square  sense  and  depolarization  Is  negligible.  Having 
obtained  ,  the  amplitude  of  the  mean  scattered  field  may  be  computed  as 


follows; 


<E  > 
s 


7^,1  jjj 


2 

(6.53) 


6.2 .3.1  Discussion  of  Re suits 

Attention  will  now  be  directed  toward  studying  the  surface  conditions 
which  permit  a  simplification  of  (6.42)  or  (6.51).  Of  particular  concern  are 
the  range  of  surface  parameters  which  result  in  a  one  dimensional  simplifica¬ 
tion  of  either  (6.42)  or  (6.51),  In  expanded  form  (6.42)  takes  the  following 
form; 


.91_ 

9  Ax 


<7;;3 

■  {-  =x  '‘2X  -  “xxB2x  <Bo-Bl>«x  -  “xyB2,)  *  i  fk  {'  S'  ''2y  '  %yB2y 

X 

■  BxyB2x}  ^  cxp|xj+J  I  dA^ 


-iBid62,d6^y 


•  2“<,K-J<^x'’x  '‘2x/^-J  S'*/''2y/']*"'r'‘lxW‘2x>«><''2y> 
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.  -J 


where 


*i-“(o^-R)(B  -kJ-(S^+R  )B„  k_  -  (S^+R  ){i-  k„  -  (k,-B  )(R  B,  +R  3,  ) 

1  '  '  0  1  '  o  1'  '  X  XX  2x  2x  '  y  yy  2y  2y  '1  o  x^2x  y  2y' 

+  (B  -B, ) (R  k-  +R  k„  )  -  R  (Bo  K  +Bo  K  ) 
o  x  2x  y  2y  xy  2y  2x  2x  2y 

When  (6.54)  can  be  essentially  reduced  to  an  integration  over  the  height 
transform  variable  (B^)  >  it  la  reasonable  to  expect  that  the  equivalent  planar 
current  is  approximately  independent  of  the  surface  slopes.  If  this  is  true 
then  Jq^^l»^2x’^2y^  aeaumes  the  following  approximate  form; 


(6.55) 


arid  substituting  this  simplification  in  (6.54)  yields 


<V2x*V2y>}  “'t  ''®oj  '*®1 

'  2H„[C„  -  ]  kj^^a  -  J  Cy  Sjf  kj^/  ‘J^'kl-l'i.kKkzxWkyy) 


(6.56) 


One  way  of  satisfying  (6.56)  is  to  take  the  exponent  which  depends  upon  R^ 
and  R  to  be  very  small.  The  maximum  amplitudes  of  R  and  R  are  of  order 

y  X  y 

S  a  end  S  O  ,  respectively.  Introducing  the  normalized  height  transform  vari- 
X  y 

ables  tIq  ®  ^o^^iz’  ^1  *  ^l^^iz  *^1  "  exponent  in  question  can 

be  Ignored  provided 


k .  as  «  1 

Iz  X 


k_,  OS  «  1 
iz  y 


(6.57) 
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AHsiimiiiK  that  (6,.*)?)  is  sntisfiod,  (6.S6)  reduces  to  the  following; 


i/ 


*//[■ 


^x9Ax  ^  y'^2y9Ay  h 


] 


•  exp  (a^-RXe^-Bj^Xe^-kj)  I  exp(J  dAr^  dB^  'dBj^ 

-  2H  fc  -JC  S^k-  /2  j  C  S^k-  /2l(l)(k,-k,  )  (6.58) 

oLo  XX  2x/  y  y  2y/  1  iz' 

Equating  like  powers  of  and  k2y>  (6.58)  yields  three  equations  for  jq(3]_)  • 

Since  it  is  not  at  all  clear  that  these  three  equations  will  yield  the  same 

2  2 

jq(Bl)  ,  it  is  further  necessary  to  assume  that  <<  1  and  <<  1  . 

Thus,  (6.58)  finally  becomes 

2ir-j 


•  dBj^  «  2H^  ct  (k^-k^^)  (6.59) 

which  is  the  desired  result.  It  should  be  noted  that  although  the  solution 
of  (6.59)  will  not  satisfy  (6.58)  exactly,  the  stipulation  of  small  mean 
square  slopes  insures  that  the  difference  will  be  small. 

Tlie  approach  presented  above  for  reducing  the  dimensionality  of  (6.42) 
or  (6.51)  is  by  no  means  unique,  i.e.  there  may  well  be  many  other  equally 
valid  techniques.  However,  the  end  result  and  the  conditions  are  expected  to 
be  the  same.  The  condition  of  small  mean  square  slopes  permits  Ignoring  the 
terms  in  (6.58)  which  exhibit  a  linear  dependence  upon  k2^  and  k2y  .  Condi¬ 
tion  (6.57)  or  the  assumption  that  the  product  of  the  rms  slopes  and  the  Ray¬ 
leigh  parameter  (l^^^®)  is  very  small  is  necessary  in  order  to  eliminate  the 
coupling  between  the  height  and  slope  resulting  from  the  integral  term  in  (6.5). 
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Of  these  two  sets  of  conditions,  it  is  obvious  that  (6.5?)  is  the  more  dif¬ 
ficult  to  satisfy  because  as  increases  the  rms  slopes  (S^  and  S^)  will 

have  to  necessarily  decrease.  This  is  a  very  Important  condition  because  it 
shows  that  (6.59)  is  valid  only  when  an  increase  in  roughness  height  Is  accom¬ 
panied  by  an  increase  in  the  horizontal  scale  of  the  rouglmess.  Thus,  the 
validity  of  (6.59)  is  intimately  tied  to  the  Interplay  between  the  rms  surface 
slope  and  the  electrical  height  of  the  surface  rouglmess.  To  the  author's 
knowledge,  this  is  a  new  result  and  it  further  Illustrates  the  power  of  the 
multivariate  approach.  That  is,  reducing  the  dimensionality  of  the  multivari¬ 
ate  Integral  equation  can  always  be  directly  linked  to  ignoring,  due  to  small¬ 
ness,  some  physical  property  of  the  surface.  This  is  one  clear  advantage  that 
the  multivariate  approach  enjoys  over  the  conventional  multiple  scattering 
approach  because  it  is  most  difficult  to  translate  closure  of  the  multiple 
scattering  equations  into  an  equivalent  surface  assumption. 

Having  demonstrated  how  one  goes  about  determining  when  the  slope  effects 
in  (6,42)  or  (6.51)  can  be  Ignored,  it  is  now  possible  to  back  up  and  estimate 
more  definitively  when  the  curvature,  rate  of  change  of  curvature,  etc.,  effect 
can  be  Ignored.  This  is  a  somewhat  tedious  task  and  will  be  saved  for  future 
studies.  However,  it  is  clearly  obvious  that  the  dominant  smallness  assump¬ 
tion  is  going  to  result  from  the  statistical  coupling  or  correlation  between 
the  surface  height  and  the  order  of  surface  height  derivative  in  question. 

Regularizing  the  Ar ^-integrand  of  (6.59)  at  infinity  leads  to  the  slope 
Independent  analog  of  (6.51),  l.e. 
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I 


-  aitp|-o^(6^-Bj)(3„-kj)}]  exp  (Jk^-A^|.)<lA?|.dB^  -  rCBjA^.V®' [ '‘^1 


=  2H„C^t(kj-k,^) 


(6.60) 


where  ?(•)  is  given  by  (6.49)  and  (6.50).  In  the  limit  of  a  randomly  ele- 
2 

vated  plane  R-^a  and  (6.60)  reduces,  as  it  should,  to  the  physical  optics 
result  for  the  mean  scattered  field.  Introducing  the  normalized  variables 


=  3^  /  K 
o  o/ 


iz 


"i  '  h/'^±z 


"i  “  h/' 


'iz 


(6.60)  becomes 

k. 


iz 

2ii 


-  exp  I  -  (no-Kj^)|Jexp(jit^*Arj.)dArj.dr,^  -  r(nj^><j^,k2=0,^2“^)/  (ZlJ)' 


*  dn,  = 


(6.61) 


where  p(Ax,Ay)  =  R(Ax,Ay)/o  is  the  normalized  surface  height  correlation 

function.  If  k,  0  <<  1  then  (6.61)  may  be  approximated  as  follows; 
iz 


p(Ax,Ay)  (Hq-Hj^)  (nQ-Kj^)exp(jk^*Ar^)dAr^dn^^  -  r(rij^,K:^,k2=0,32=0)/ (2it)' 


=  2KoC^cJ,(k^^[K^-l)) 


dn. 


(6.62) 


102 


9ivf!smasii['Ziri'.-j^4 


where  the  exponential  term  In  (6.61)  containing  p(Ax,Ay)  has  been  approxi¬ 
mated  by  a  two  term  power  series  because  ^  *  Consistent  with  the 

previous  assumptions  of  small  mean  square  slope,  curvature,  rate  of  change 
of  curvature,  etc.  p(Ax,Ay)  may  be  conslaered  to  vary  much  more  slowly  with 
Ax  and  Ay  than  and  so  (6.62)  becomes 


J (V'l)  V"l>  ic 

•  exp  (k^^O)^(ri^-n^)(n^-Kj^)|dn^  -  r(n^,Kj^,k2=0,B2=0)/(2lT)^^  dTij^ 


(6.63) 


Since  the  integral  term  is  multiplied  by  ,  it  is  much  smaller  than  r(») 

because  k^^O  «  1  ;  thus,  for  small  Rayleigh  parameter  the  integral  equation 


for  1  (k.  n, )  becomes 
''q  iz  1 


»  2H  C  't>(k  [K.-l]) 

o  o  iz  X 


(6.64a) 


or  in  terms  of  and  kj^ 


^  J “  r(Bi.kj,k2*°  *^^1 


»  2H^  ^ (k^-k^^)  (6.64b) 

i’liei'o  are  some  iitLeresi  iiig  eonHe<(uei(eeH  ul  ((i.(»4)  wlih'li  will  now  In*  i'oie 
sidered.  First,  one  way  of  achieving  a  small  Rayleigh  parameter  is  to  let  0^ 
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Pi  : 
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f,  ;' 

r>: 


t  i 


n 

u  i 


f,  ■ ; 
h  * 
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be  very  near  grazing  Incidence  (0^^  ■  Tr/2)  .  Near  grazing  incidence  is  where 

approximate  optical  theories  predict  a  strong  dependence  on  the  mean  square 

slope  due  to  self-shadowing  [^,5].  However,  (6.64)  exhibits  no  dependence 

upon  the  surface  slopes  and,  in  fact,  is  determined  entirely  by  the  Rayleigh 

parameter  ‘  This  discrepancy  between  (6.64)  and  approximate  optical 

theories  does  not  appear  to  be  a  consequence  of  any  of  the  simplifications 

leading  to  (6.64)  because  these  same  simplifications  are  inherent  in  the 

approximate  theory  (except  for  k.  OS  «1  which  is  trivially  satisfied 

xz 

near  grazing  incidence).  The  lack  of  dependence  of  (6.64)  upon  the  surface 
slope  variances  is  also  in  agreement  with  rigorous  boundary  perturbation 
theory  which,  in  turn,  is  accurate  for  k  o  <<  1  .  Thus,  (6.64)  clearly 

L  Z 

establishes  the  inaccuracy  of  shadow  corrected  optical  approximations  for 
near  grazing  incidence  or,  more  correctly,  for  small  Rayleigh  parameters,  at 
least  for  the  case  of  coherent  scattering.  One  possible  reason  for  the  fail¬ 
ure  of  the  shadow  corrected  optical  approximation  is  that  it  assumes  that  the 
current  Induced  on  the  surface  is  zero  on  the  shadowed  parts  of  the  surface. 
However,  this  is  at  complete  variance  with  (6.64)  which,  for  k.  O  suffici- 

XZ 

ently  sinall,  predicts  no  shadowing  of  the  incident  field.  In  fact  (6.64) 

shows  that  for  a  sufficiently  small  Rayleigh  parameter,  the  average  scattered 

field  appears  to  result  from  a  randomly  elevated  plane. 

The  average  scattered  field  for  the  case  of  k.  0  <<  1  can  be  obtained 

Iz 

by  first  solving  (6.64)  for  then  convolving  this  result  with 

ct>(k  -3.),  k  =  -k,  ,  to  find  the  average  scattered  field.  There  is,  how- 

S  Z  X  S  Z  1-  z 

ever,  a  more  direct  approach  to  computing  <E  >  which  results  from  the  special 

s 

form  of  (6.64).  Writing  !’(•)  in  its  Integral  form  as  given  by  (6.46)  yields 
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(6.65 


Since 


2  2 

4.(kj-Bj)  (B„-kj)].  «xp|-  |-(kj-B„)^  -  r(ei-6o>^} 


-  4>(kj-B„)«B„-Bi) 


="s»„)>  -  ^  J l,(BjWB„-B^)  dBj 


(6.65)  may  be  rewritten  as  follows 


<E^(k^)>.2H^C„<Kkj-k^^)+  j'<B,(B„»4>(kj-B„)l^(£,  .B„)dB„/:2p)' 


(6.66 


Furthermore,  since 


where  e  is  a  small  positive  quantity,  (6.66)  reduces  to  the  following  sin 
guiar  Integral  equation  of  the  second  kind  for  <H  >  ; 


<E  (k-)>  »  2Ii  C 
s  i  o  o 


<l)(k-k  )+i  I 


<E  (B  )>  dB 
s  o  o 


(6.67 


'Hi is  eqiiutLon  bears  a  remarkable  resemblance  to  a  result  obtained  by  DeSanto 

2 

[8,10],  wliicb  lie  stated  was  valid  In  the  limit  of  R-^O  .  In  comparing  (6.67) 

with  UeSanto's  result,  there  are  two  factors  that  should  be  noted.  First, 

DeSanto' s  equation  has  a  in  the  denominator  of  the  integral  rather  than 

in  the  numerate  This  difference  appears  to  result  from  the  manner  In  which 

he  treats  the  singularity  in  Bg/DC  .  In  (6.67)  the  3  In  the  numerator 

o 

results  very  simply  from  the  following  transform  relationship; 

For  k.  (7  <<  1  ,  tills  difference  will  have  no  essential  effect  upon  (k  )>  . 
A  second  point  to  be  noted  is  the  fact  that  DeSanto' s  result  corresponding  to 
(6.67)  is  actually  based  upon  tlie  assumption  that 

Ar^.  -  0 

R(Arj.)  =  (6.68) 

0  A  r^.  0 

2 

rather  than  [<->0“"  as  stated  in  [9]  and  [lOj.  Since  the  correlation  func¬ 
tion  in  (6.68)  differs  from  zero  only  over  a  domain  of  zero  measure,  substi¬ 
tution  of  (6.68)  In  (6.60)  will  thus  yield  (6.67).  In  fact  (6.68)  was  dis¬ 
covered  by  searching  for  the  form  of  R(Ar^)  which  would  reduce  (6.60)  to 
(6.6/).  However,  it  must  be  remembered  that  (6.60)  is  based  upon  the  assump¬ 
tion  that  the  variances  of  V”r;»n=l  ,2  , . . .  ,  are  all  very  small  and  it  is  not 
Immediately  obvious  that  (6.68)  satisfies  these  conditions.  Thus,  although 
(6.6  7)  can  be  obtaint='d  in  the  same  manner  as  DeSanto  derived  his  correspond¬ 
ing  result,  there  remains  some  question  as  to  the  meaning  of  the  expressions 
for  •  Aside  from  these  minor  differences,  it  is  most  encouraging  that 

these  two  significantly  different  mathematical  approaches  give  rise  to  veiry 
similar  equations  for  the  average  scattered  field. 


106 


When  the  Rayleigh  parameter  becomt's  mod^.-rate  to  large,  (6.67)  can  no 
longer  I  .  rigorously  justltied  and  (6.60)  must:  be  solved,  blnce  the  only 
major  difference  between  (6.67)  and  (6.60)  la  the  appearance  in  (6.60)  of 
the  tern  involving  integrations  over  Ar^.  and  3^  »  it  would  appear  that 
(6.60)  should  be  readily  aoluable  also.  Unfortunately,  no  an 'iytlcai  approxi¬ 
mations  have  been  found  for  this  term.  Since  tliis  term  is  roughly  equivalent 
to  the  average  of  the  i;-derivative  of  the  Green’s  frnction,  it  is  extremely 
importamt  to  the  determinatiixa  of  or  •  Quite  obviously,  future  efforts 

on  this  problem  should  concentrate  on  the  approximate  analytical  evaluation 
of  this  term.  If  analytical  approaches  prove  fruitless,  then  numerical  inte¬ 
gration  techniques  should  also  be  consid  red. 

6 ■ 3  Summary 

The  purpose  of  this  section  is  to  develop  an  alternate  approach  to  the 
problem  of  coherent  scattering  from  a  perfectly  conducting  rough  surface. 

Since  the  conventional  multiple  scattering  approach  is  difficult  to  interpret 
in  terms  of  the  statistical  properties  of  the  surface,  it  is  desirable  to  have 
an  approach  which  clearly  shows  the  dependence  of  the  mean  scattered  field 
upon  the  various  surface  parameters,  ‘ilie  approach  developed  in  this  section 
leads  to  a  single  integral  equation  of  infinite  dimensions  for  a  function 
which  when  convolved  with  the  Joint  characteristic  function  for  the  height, 
slope,  curvature,  etc.  leads  to  the  coherent  field.  The  Integral  equation 
has  infinite  dimension  because  the  mean  scattered  field  depends  upon  all  order 
derivatives  of  the  surface  height.  Particular  attention  is  given  to  a  Gaussian 
surface  which  has  negligible  curvature,  rate  of  change  of  curvature,  etc.  It 
is  shown  how  one  links  reducing  the  dimensionality  of  the  integral  equation 
to  conditions  on  the  surface  parameter.  In  particular  It  is  demonstrated  that 


107 


the  equivalent  planar  current  Induced  on  the  surface  can  be  taken  to  be  Inde¬ 
pendent  of  the  surface  slopes  whenever  the  slope  variances  are  small  and  the 
product  of  the  Kaylelgh  parameter  and  the  rms  slope  is  small.  For  near  graz¬ 
ing  Incidence  it  is  found  that  there  is  no  justification  for  the  inclusion 
of  a  shadowing  function.  Comparing  these  results  with  those  obtained  from 
multiple  scattering  theory  shows  a  great  deal  of  agreement  and.  In  addition, 
some  further  Insight  Into  certain  limiting  approximations. 

tn  summary  the  multivariate  approach  introduced  here  has  the  following 
positive  qualities; 

•  it  is  very  straightforward, 

•  it  is  sufficiently  general  to  cover  any  surface  height 
distribution, 

•  it  clearly  shows  the  dependence  of  the  scattered  field  on 
measurable  surface  parameters, 

•  one  can,  in  a  straightforward  mamier,  deduce  when  the  various 
higher  order  surface  height  derivatives  are  important, 

•  tlie  technique  can  be  extended  to  tlie  general  vector  case,  and 

•  the  basic  approach  can  also  be  used  to  determine  the  incoherent 
field. 

The  only  negative  aspect  of  the  approach  is  that  it  will  require  the  evalua¬ 
tion  of  some  very  complicated  and  difficult  integrations  which,  in  essence, 
stem  from  averaging  spatial  derivatives  of  the  Green's  function. 
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